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FORWARD  -  EXECUTIVE  SUMMARY 


Development  of  rapid  maneuvering  control  algorithms  for  Army  supersonic  missiles 
at  low  altitudes  using  low  blowing  ratio  reaction  jets  is  complicated  by  nonlinear, 
counterintuitive  surface  interactions,  hot  gas  effects  and  poor  predictability  using  modem 
computation  fluid  dynamic  codes,  where  turbulence  modeling  is  a  key  factor  limiting  the 
accuracy.  The  objectives  of  this  research,  which  were  directly  tied  to  the  Army  low  blowing 
ratio  reaction  control  jet  requirements,  were  documentation  and  improved  understanding  of 
(1)  the  surface  flow  interactions  with  hot  gas  effects  and  (2)  the  turbulent  flow  structure  in 
the  plume.  An  experimental  analysis  of  the  mean  and  turbulent  flow  properties  for  sonic  and 
Mach  3  injection  into  a  Mach  3  freestream  was  accomplished.  The  range  of  conditions 
included:  MPR  =  0.005-0.09,  Pj  /  />*,=  0.6-19.0,  Tj  / Tx  =  1.0-13.6  (helium  was  used  to 

simulate  the  highest  temperature),  y  =  1.4  and  1.67,  and  peb  /  pK  =  0.19-6.15 
(overexpanded,  perfectly  expanded  and  underexpanded  jets  were  tested).  Numerical 
simulations  were  not  accomplished  in  the  present  effort.  However,  a  full  characterization  of 
the  high-speed  jet  interaction  flowfield  over  a  blended  body  missile  fuselage  (sharp-coned 
cylinder)  was  experimentally  achieved,  where  specific  pragmatic  and  scientific  flowfield 
conclusions  were  drawn.  For  example,  injection  temperature,  for  a  given  overexpanded 
MPR,  had  a  small  but  discernible  effect  on  the  boundary  layer  separation  distance.  The 
influence  of  the  jet  plume  on  the  surface  pressure  distribution  in  the  wake  region  of  the  flow 
was  found  to  depend  on  MPR,  where  the  wake  region  developed  further  downstream  in  terms 
of  distance  nondimensionalized  by  the  jet  exit  port  diameter.  The  effective  backpressure 
concept  was  an  accurate  indicator  of  the  jet  flow  condition  (overexpanded,  underexpanded  or 
perfectly  expanded).  The  upstream  peak  pressures  associated  with  lambda  shock  and  the 
horseshoe  vortex  increased  with  increasing  MPR.  The  wake  pressures  decreased  with 
increasing  MPR.  The  interaction  force  increased  almost  linearly  with  increasing  MPR.  The 
amplification  factor  was  inversely  proportional  to  MPR.  In  general,  progressing  downstream, 
the  turbulence  levels  formed  three  distinct  regions  across  the  plume.  Region  1  was  defined  as 
the  upper  region  of  the  plume,  where  the  turbulence  levels  were  moderate.  Region  2  was 
defined  as  the  high  turbulent  level  region  located  along  the  lower  portion  of  the  plume  near 
the  missile  fuselage.  Region  3  was  define  as  the  relatively  low  turbulence  region  along  the 
plume  core,  where  in  the  upper  half,  the  turbulence  levels  were  very  low,  and  moving  down 
towards  the  missile  fuselage,  the  levels  increased  to  those  in  region  2.  The  present  data 
suggest  a  strong  qualitative  correlation  between  the  measured  shear  strain  rate  field  and  the 
turbulence  across  the  plume,  which  supports  the  usage  of  the  Boussinesq  turbulence  model 
approximation. 
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1.  PROBLEM  STATEMENT  AND  SCIENTIFIC  PERSONNEL 

1.1  Problem  Statement 

1.1.1  Army  Motivation 

Jets  have  the  potential  to  be  an  effective  means  of  controlling  aerodynamic 
characteristics  of  high-speed  vehicles  and  weapon  systems.  One  Army  application  is  the 
use  of  reaction  control  jets  for  rapid  maneuvering  of  high-speed  intercept  missiles. 
Hypersonic  interceptor  designs  typically  include  several  lateral  jets  used  for  divert  and/or 
altitude  control.  Air-to-air  interceptors  usually  require  a  continuous  and  long  duration 
reaction  control  jet  for  accuracy  and  rapid  maneuverability.  In  addition,  confined  space 
firing  (e.g.,  anti-tank  or  Guerilla  warfare)  missile  control  systems  are  required  to  be 
effective  at  low-speed  with  good  maneuverability  and  very  short  response  time 
[Champigny  and  Lacau  (1994)]. 

Accurate  aerodynamic  characterization  of  the  jet  interaction  with  the  missile  body 
is  essential  for  the  design  of  an  autopilot.  One  important  parameter  of  jet  interaction  is  the 
amplification  factor,  which  is  defined  as  the  ratio  of  the  effective  force  to  the  vacuum 
thrust.  Vacuum  thrust  is  the  net  reaction  force  created  when  a  jet  is  fired  in  a  vacuum;  this 
force  is  equivalent  to  the  thrust  of  a  rocket  engine.  However,  in  a  sensible  atmosphere, 
the  net  force  produced  between  the  jet  and  free  stream  can  be  higher  or  lower  than  the 
vacuum  thrust  of  the  rocket  engine.  In  the  extreme  case,  the  net  force  on  an  axisymmetric 
body  can  be  the  opposite  of  what  is  expected,  a  situation  known  as  control  reversal 
[Chamberlain  et  al  (2000)].  When  a  jet  is  first  fired  into  a  supersonic  crossflow,  a  shock 
wave  upstream  of  the  jet  is  created.  The  high  pressures  associated  with  the  shock  causes  a 
positive  amplification.  A  portion  of  the  shock  that  wraps  around  the  missile  fuselage 
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creates  a  high-pressure  region  under  the  axisymmetric  body  causing  a  de-amplification 
on  the  net  force.  Also,  the  low-pressure  wake  region  downstream  of  the  jet  causes  de¬ 
amplification  factor  is  further  decreased.  If  the  total  de-amplification  is  higher  than  the 
total  amplification,  the  net  force  will  be  opposite  than  what  is  anticipated.  The  most 
severe  de-amplification  happens  at  low  altitudes,  low  Mach  numbers,  and  low  angle  of 
attack. 

At  high  altitude,  typical  reaction  control  jets  produce  several  thousand  pounds  of 
thrust  in  less  than  10  ms  [Roger  (1999)].  The  forces  involved  at  high  altitude  are 
fundamentally  of  action-reaction  type  in  the  usual  Newtonian  sense.  However,  at  lower 
altitudes,  jet  interaction  is  confined  to  be  near  the  surface,  and  hence  produces  very 
complicated  nonlinear  surface  interactions,  which  can  adversely  affect  the  vehicle 
control.  These  nonlinear,  and  often  unsteady,  interactions  must  be  understood  in  order  to 
develop  effective  and  accurate  vehicle  control  schemes. 

A  number  of  parameters  have  been  listed  that  can  affect  thruster  control 
effectiveness.  A  subset  of  salient  parameters  are  [Roger  (1999)]:  the  freestream  Mach 
number  (M*,),  the  jet  Mach  number  (Mj),  the  jet  Momentum-Parameter-Ratio  (MPR),  the 
jet  pressure  ratio  (Pj/Poo),  the  jet  temperature  ratio  (T/T*,),  the  Reynolds  number  (Re),  and 
the  jet  gas  ratio  of  specific  heats  (yj).  In  addition  to  these  parameters,  jet  geometry,  and 
chemical  reactions  within  the  plume  were  also  identified  as  important,  especially  for  flow 
separation. 
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1. 1.2  Army  Research  Challenges 

Flight  control  of  supersonic  Army  missiles  at  low  altitudes  using  low  blowing 
ratio  reaction  jets  is  complicated  due  to  (1)  nonlinear,  counterintuitive  surface 
interactions,  (2)  hot  gas  effects  and  (2)  poor  predictability  using  modem  computation 
fluid  dynamic  codes,  where  turbulence  modeling  is  a  key  factor  limiting  the  accuracy 
[Doligalski  (1997)].  Furthermore,  it  has  been  documented  [Champigny  and  Lacau  (1994), 
Chamberlain  et  al  (1999),  Dash  et  al  (1999),  Graham  et  al  (1999),  Hsieh  (1999),  Thoenes 
(1999),  and  Roger  (1999)]  that  jet  interaction  experimental  data  is  lacking,  and  high 
quality  systematic  experimental  studies  to  build  a  database  are  needed.  Numerous  three- 
dimensional  vortex  features,  unsteadiness,  turbulence,  large  thermal  gradients,  and  real 
gas  effects  characterize  the  flowfield  associated  with  a  jet  into  a  high-speed  crossflow. 
These  features  make  the  flowfield  very  complicated.  A  detailed  description  of  the  up-to- 
date  research  contributions  on  the  flowfield  is  provided  in  Chapter  2.  The  poor  flowfield 
understanding  and  predictability  are  underscored  by  reported  computational  fluid 
dynamic  calculations  that  predicted  the  normal  force  in  the  wrong  direction  [Doligalski 
(1997)].  Appropriate  turbulence  modeling  is  the  key  factor  that  limits  numerical 
accuracy. 

Currently,  engineers  and  scientists  rely  on  an  approximate  form  of  the  governing 
equations  of  motion  to  predict  the  overall  meanflow  character  of  these  flows.  The 
Reynolds  and  Favre  averaged  forms  of  the  Navier-Stokes  equations  are  universally  used 
to  compute  low-speed  and  high-speed  flows,  respectively.  Some  recent  studies  [Graham 
and  Weinacht  (1999),  Hsieh  (1999),  Kikumoto  and  Sentoh  (2001),  Graham  and  Weinacht 
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(2000)]  (using  Reynolds  averaged  N-S  technique,  Baldwin-Lomax,  k-s,  GASP,  PARCH, 
etc)  showed  that  existing  turbulence  models  fail  near  the  recirculation  region  in  front  of 
the  jet  and  near  the  recompression  shock.  To  correct  these  shortcomings,  there  is  a  need 
for  high  quality,  high  fidelity  experimental  data  that  will  allow  for  verification  of  current 
and  future  computational  fluid  dynamic  results  [Margason  (1993)]. 

1.1.3  Research  Objectives 

The  overarching  objective  of  the  present  research  is  to  improve  the  understanding 
and  predictability  of  high-speed  high-temperature  jet  interaction  flowfields.  Specific 
Army  requirements  are  low  Mach,  low  altitude,  and  low  blowing  ratio  operating 
conditions  missile  flight  control.  Therefore,  the  following  objectives  were  established  for 
this  research. 

1.  Documentation  and  improved  understanding  of  the  missile  surface 
interactions  for  low  MPR  jets,  with  elevated  injection  temperature. 

2.  Documentation  and  improved  understanding  of  turbulent  flow  structure  of  the 
jet  plume  for  the  low  MPR  jet  interaction  flows. 

1. 1.4  Research  Approach 

A  three-year  research  program  was  fully  dedicated  to  experimental 
characterization  of  the  jet  interaction  flowfield  over  a  3-D  blended  missile  fuselage 
model.  A  full  gamut  of  experimental  techniques  was  used  to  characterize  the  mean  and 
turbulent  flow  properties.  The  tunnel  was  operated  to  simulate  Mach  3.0  flight  at  an 
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altitude  of  17-km.  The  range  of  jet  parameters  studied  here  is  listed  in  Tables  1,  and  the 
freestream  flow  conditions  are  given  in  Table  2. 

In  summary,  surface  oil  flow  visualization  and  pressure  sensitive  paint 
experiments  were  performed  to  provide  a  global  picture  of  the  flow/surface  interaction 
and  aerodynamic  loading.  Schlieren/shadowgraph  techniques  were  used  to  quantify  the 
shock  structure.  Instantaneous  and  time  averaged  Mie-scattering  imaging  was  used 
characterize  the  jet  plume  trajectories.  Mean  flow  mapping  was  accomplished  using  a  5- 
hole  pressure  probe  and  particle  image  velocimetry.  Particle  image  velocimetry  also 
provided  documentation  of  the  turbulent  flow  structure  in  the  plume. 

1.2  Scientific  Personnel 

The  PI  for  this  project  was  Dr.  Rodney  D.  W.  Bowersox,  Associate  Professor  of 
Aerospace  Engineering.  This  project  supported  the  PhD  dissertation  research  of  Mr. 
Zakaria  Mahmud,  who  successfully  defended  his  dissertation,  in  the  Department  of 
Aerospace  Engineering  and  Mechanics  at  the  University  of  Alabama,  on  November  8, 
2002. 
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2.  BACKGROUND  REVIEW 


2.1  Continuous  Injection  into  a  Low-Speed  Crossflow 

The  flowfield  for  injection  in  a  cross  flow  research  has  received  considerable 
attention  because  of  the  many  important  practical  applications  in  both  military  and  civil 
disciplines.  Margason  (1993)  summarized  the  evolution  of  jet  into  a  cross  flow  research 
for  the  time  span  between  1932  and  1993.  Over  300  articles  were  included  in  this  review. 
A  brief  synopsis  of  the  salient  features  is  presented  here  with  the  inclusion  of  results  from 
recent  work. 

The  numerous  studies  on  transverse  jets  into  low-speed  crossflows  [Keffer  et  al 
(1963),  Kamotani  (1972),  Moussa  (1977),  Perry  (1993),  Fric  and  Roshko  (1994),  Findlay 
(1999)]  provide  detailed  documentation  of  the  mean  flowfield  structure  and  qualitative 
insight  into  the  vortical  flow  features.  The  flowfield  structure  for  injection  into  a  low- 
speed  crossflow  is  shown  schematically  in  Fig.  1.  There  appears  to  be  at  least  four 
vortical  flow  motions  created  by  the  transverse  jet  in  a  crossflow  [Fric  and  Roshko 
(1994)].  As  the  jet  emerges  into  the  crossflow,  it  is  curved  downstream  by  the  crossflow. 
The  jet  plume  acts  as  a  cylindrical  obstacle  to  the  crossflow  and  creates  an  adverse 
pressure  gradient,  which  eventually  generates  a  wraparound  horseshoe  vortex  on  the 
surface.  Then,  the  jet  plume  evolves  into  a  “kidney-bean”  shape,  which  is  due  to  the 
counter  rotating  vortex  pair.  This  counter  rotating  vortex  pair  dominates  the  downstream 
flow  [Kamotani  (1972)].  The  origin  of  this  vortex  pair  is  believed  to  be  near  the  jet  exit 
and  caused  by  the  lateral  shearing  along  the  plume  edges  and  the  pressure  gradient  due  to 
turning  the  flow  into  the  freestream  direction.  Another  ring-like  vortex  structure  appears 
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in  the  wake  region.  These  alternating  vertical  vortices,  in  many  papers  named  as  the  wake 
vortices,  give  the  appearance  of  a  “Karman-vortex-street” [Perry  et  al  (1993)]. 

Yasuhiro  and  Isaac  [Kamotani  et  al  (1972)]  used  hot-wire  anemometry  to  map  the 
velocity,  temperature,  and  turbulence  intensities  for  a  low-speed  JI  flow.  They  discerned 
that  the  momentum  ratio  significantly  affects  the  jet  velocity,  and  turbulence  intensity; 
while  the  entrainment  process  is  independently  controlled  by  the  normal  and  parallel 
components  of  the  velocity.  Perry  et  al  [Perry  et  al  (1993)]  studied  the  entrainment  and 
mixing  processes  of  a  low-speed  jet  in  a  cross  flow.  The  results  from  the  flying  hot-wire 
and  flow  visualization  revealed  that  there  exists  many  vortex  systems  and  they  are  inter¬ 
connected;  the  karman  vortex  street  like  structure  come  from  the  vorticity  generated  at 
the  wall.  Findlay  et  al  [Findlay  et  al  (1999)]  employed  three-component  LDV  to  study  the 
flowfield  characteristics.  The  paper  further  proved  the  accepted  flowfield  characteristics. 

2.2  Continuous  Injection  into  a  High-Speed  Crossflow 

The  compressibility  effects  associated  with  injection  into  a  high-speed  crossflow 
create  additional  flowfield  complications.  An  extensive  literature  review  for  high-speed 
jet  interaction  flow  was  performed  as  part  of  this  study.  Over  100  articles  were  examined. 
A  synopsis  of  the  flowfield,  built  from  this  review  for  flat  and  axisymmetric  surfaces,  is 
described  below. 

2.2.1  Jets  Issuing  from  a  flat  Surface  into  a  Supersonic  Stream 

Many  of  the  mean  flow  features  of  injection  into  supersonic  crossflows  are  similar 
to  those  described  above  for  a  low-speed  crossflow.  However,  compressibility  creates 
additional  features  that  do  not  have  incompressible  counterparts  [Schetz  and  Billig 
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(1966),  Schetz  et  al  (1967),  Billig  et  al  (1971),  Heister  and  Karagozian  (1989),  Grasso 
and  Magi  (1995),  Santiago  et  al  (1995)  and  Chenault  et  al  (1998)].  Figure  2(a)  presents 
the  accepted  3-D  structure  of  a  highly  under  expanded  jet  into  a  supersonic  cross  flow  on 
a  flat  surface.  The  transverse  jet  causes  an  obstruction  to  the  freestream  flow  and  a  bow 
shock  is  created  upstream  of  the  jet;  typically  termed  as  the  interaction  shock.  For 
injector  configurations  where  the  injector  diameter  is  greater  than  the  boundary  layer 
thickness,  a  separation  bubble  and  a  lambda  shock  form  upstream  of  the  injector  port 
[Schetz  and  Billig  (1966)].  For  most  applications,  the  jet  is  operated  in  an  underexpanded 
condition.  Hence,  the  jet  fluid  undergoes  a  Prandlt-Meyer  expansion,  recompresses 
through  a  barrel  shock,  and  then  terminates  into  a  normal  shock,  which  is  called  the  Mach 
disk.  As  the  jet  plume  curves  downstream  due  to  the  crossflow,  the  shape  of  the  plume  is 
transformed  into  the  kidney  bean  shape,  which  is  attributed  to  the  counter  rotating  vortex 
pair.  As  the  freestream  flow  wraps  around  the  jet,  a  recompression  shock  is  created. 
Chenault  et  al  (1999)  suggested  that  the  recompression  shock  was  a  new  mechanism 
adding  to  the  strength  of  the  vortex  pair. 

The  effective  backpressure  plays  an  important  role  on  the  structure  of  the  jets  in  a 
supersonic  crossflow  [Schetz  et  al  (1967)].  In  addition  to  underexpanded  jets,  the  present 
research  also  examined  overexpanded  and  perfectly  expanded  cases.  For  the 
overexpanded  case,  the  flowfield  structure  is  slightly  different.  Here,  the  jet  plume 
terminates  through  an  oblique  shock  cone.  The  angles  of  the  conical  shock  system  depend 
on  the  blowing  ratios.  Heister  and  Karagozian  (1989)  studied  the  perfectly  expanded 
transverse  jets  into  the  supersonic  crossflow,  where  terminating  shocks  were  absent. 
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In  1960s,  the  flowfield  structure  was  documented  using  schlieren  photographs, 
mean  flow  measurements,  concentration  measurements,  and  mass  flow  rates.  Schetz  et  al 
[Schetz  et  al  (1967),  Schetz  et  al  (1968)]  investigated  the  relationship  between  the  highly 
under  expanded  jet  into  a  supersonic  (Mach  2.1)  freestream,  as  well  as  for  the  jet  issuing 
into  a  quiescent  medium  to  determine  the  penetration  height  and  the  horizontal 
displacement  of  the  mach  disk.  The  papers  concluded  that  for  a  fixed  Mj  and  fixed  jet 
geometry,  the  displacement  of  the  Mach  disk  varies  as  Pj1/2;  and  the  penetration  height  is 
proportional  to  the  half  power  of  the  mass  flow  rate  at  each  Mach  number.  The 
penetration  height  increased  with  increasing  jet  Mach  number,  however  the  effect  is  not 
very  significant  beyond  Mach  2.0.  Billig  et  al  (1971,  1994)  later  developed  a  model  to 
account  for  different  jet  injection  angles  and  the  turbulent  mixing  after  the  Mach  disk 
using  a  modified  effective  backpressure  concept. 

The  effects  of  compressibility  on  free  shear  layers  are  also  important  area  of 
research  [Samimy  et  al  (1990),  Samimy  et  al  (1992),  Papamoschou  and  Roshko  (1988), 
Bradshaw  (1966),  Papamoschou  and  Roshko  (1986),  Samimy  et  al  (1989)].  The  concept 
of  the  convective  Mach  number  is  well  agreed  upon  to  be  the  most  effective 
compressibility  factor.  In  summary,  these  papers  documented  the  non-dimensionalized 
shear  layer,  vorticity  and  momentum  growth  rates  as  a  function  of  convective  Mach 
number,  ranging  from  subsonic  to  supersonic. 

McCann  and  Bowersox  (1996)  investigated  air  injection  (Mach  1.8)  at  an  angle  of 
25  degree  into  a  freestream  of  Mach  2.9.  Using  cross  film  anemometry  and  conventional 
mean  flow  probes,  the  paper  provided  detailed  compressible  turbulence  and  mean  flow 
measurements.  The  turbulent  kinetic  energy  had  two  peaks,  which  were  nearly  co-located 
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with  the  counter  rotating  vortex  pair.  The  compressibility  terms  ( up'v'/puV  and 
up'w'/ pu'w')  were  also  found  to  be  significant,  and  were  estimated  as  67.0  -  75.0%  of 
the  total  Reynolds  shear  stress.  Bowersox  (1996)  further  investigated  the  anisotropy 

coefficients  (e.g.;  A  =  ^j(pv)'2  ,  Ax_  and  Azv )  across  plume  for  the  same 

parametric  conditions.  The  distribution  of  A_v  was  asymmetric,  with  the  maximum  in  one 

vortex  and  the  minimum  in  the  other.  The  results  of  the  experiment  were  later  simulated 
with  the  second  order  Reynolds  stress  turbulence  modeling  [Chenault  et  al  (1999)].  The 
analysis  was  physically  consistent,  where  accurate  prediction  for  both  the  mean  and 
turbulent  flow  quantities  was  achieved. 

Santiago  et  al  (1995)  mapped  the  instantaneous  velocity  distribution  of  sonic 
transverse  jet  into  a  Mach  1.6  free  stream.  A  two-component  frequency  pre-shifted  LDV 
system  was  used  to  provide  over  2200  measurement  locations  in  the  transverse  midline 
plane.  A  potentially  important  suggestion  from  that  study  was  that  for  underexpanded 
jets,  counter  sinking  the  jet  exit  holes  to  let  it  fully  expand  before  it  encounters  the  cross 
flow  might  increase  the  penetration  height  and  mixing. 

Gruber  et  al  (1997)  in  1997  used  pressure  sensitive  paints  to  measure  surface 
pressure  in  a  supersonic  transverse  injection  flowfields.  Both  circular  and  elliptical 
nozzles  were  investigated  at  different  momentum  flux  ratios  (0.99  to  3.95),  and  it  was 
found  that  momentum  flux  ratios  as  well  as  injector  geometries  have  significant  effects 
on  the  surface  pressure  distribution.  The  research  was  conducted  with  two  different 
pressure  sensitive  paints,  where  the  PtOEP  based  paint  was  found  to  be  more  temperature 
sensitive  than  the  pyrene  based  paint. 
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2.2.2  Jets  Issuing  from  a  Body  of  Revolution  into  a  Supersonic  Stream 

The  flowfield  structure  of  a  jet  issuing  from  a  curved  body  of  revolution  is 
qualitatively  similar  to  that  over  a  flat  surface.  The  quantitative  information  varies  due  to 
the  three-dimensional  effect  of  the  curved  surface.  Figure  2(b)  presents  the  3-D  structure 
of  a  jet  emerging  from  an  axisymmetric  body  (this  sketch  was  based  on  the  present 
study). 

A  number  of  studies  [Champigny  and  Lacau  (1994),  Chamberlain  et  al  (1999), 
Dash  et  al  (1999),  Graham  and  Weinacht  (1999),  Hsieh  (1999),  Thoenes  (1999),  Roger 
(1999),  Naumann  et  al  (1993),  Kikumoto  et  al  (2001),  Graham  and  Weinacht  (2000), 
Brandeis  and  Gill  (1998),  Srivastava  (1997),  Srivastava  (1998)  and  Srivastava  (1999)] 
investigated  the  influence  of  jets  issuing  from  a  curved  surface  into  a  supersonic  cross 
flow.  Champigny  and  Lacau  [(1994)]  reported  on  the  importance  of  lateral  jet  control  for 
tactical  missiles  and  the  associated  difficulties  with  the  wind  tunnel  testing.  The  hot  gas 
effects  were  of  fundamental  concern  for  a  wind  tunnel  test  in  comparison  to  the  reactant 
jet  fluids  at  flight  conditions.  It  was  stated  that  low  molecular  weight  gas  (e.g.;  helium) 
could  improve  the  simulation,  where  for  the  same  jet  momentum  parameter  ratio  (MPR), 
the  jet  velocity  is  much  higher  than  the  cold  air  jet.  As  a  result,  helium  at  ambient 
temperatures  can  be  used  to  simulate  air  at  approximately  2000K. 

Graham  and  Weinacht  (1999)  numerically  simulated  lateral  control  jets  in 
supersonic  flow  for  an  axisymmetric  body.  Surface  pressure  measurements  downstream 
of  the  jet  were  found  to  be  in  good  agreement  with  the  experimental  results.  Using  the 
existing  Reynolds  averaged  Navier-Stokes  computational  technique,  the  paper  also 
concluded  that  the  flight  velocity  and  jet  mass  flow  appeared  to  be  the  most  significant  on 
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the  force  and  moments  control,  while  nose  shapes  and  angle  of  attack  produced  relatively 
small  effects.  In  a  second  numerical  study,  Hsieh  (1999)  found  that  the  thrust  coefficient 
and  jet  Mach  number  should  be  identical  in  order  to  have  the  satisfactory  numerical 
simulation  of  missile  control  jets. 

Roger  (1999)  summarized  nearly  fifty  years  of  research  on  flight  and  wind  tunnel 
tests  of  thruster  control  for  supersonic  and  hypersonic  interceptors.  Some  of  the  salient 
findings  were:  (1)  for  90°  jets,  there  is  essentially  no  change  in  amplification  factor  with 
jet  Mach  number;  (2)  attitude  moment  amplifications  factors  are  always  positive  and  can 
be  larger  than  three,  while  lateral  amplification  factors  can  be  negative;  (3)  force  and 
moment  measurements  can  be  affected  if  the  jet  plume  interacts  with  the  wind  tunnel  wall 
boundary  layer,  and  (4)  a  methodology  to  determine  jet  interaction  performance  for 
modem  hit-to-kill  designs  is  lacking. 

Recently,  it  was  reported  that  Reynolds  averaged  CFD  surface  pressure 
distribution  predictions  were  within  the  experimental  accuracy  [Graham  and  Weinacht 
(2000)  and  Kikumoto  and  Sentoh  (2001)].  However,  both  of  these  studies  showed 
discrepancies  in  the  recirculation  region  in  front  of  the  jet  and  at  the  end  of  the  body  near 
the  recompression  shock. 
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3.  FACILITIES  AND  INSTRUMENTATION 


3.1  Wind  Tunnel 

This  research  project  started  at  a  very  preliminary  stage,  where  as  a  first  step  the 
wind  tunnel  was  designed  and  constructed  in-house.  A  schematic  drawing  of  the  wind 
tunnel  is  shown  in  Fig.  3(a).  Figure  3(b)  shows  a  photograph  of  the  wind  tunnel.  The  free 
stream  flow  conditions  are  listed  in  Table  2.  A  detailed  specification  of  the  compressors, 
storage  tank,  nozzle  block,  models,  test  section,  and  diffusers  are  provided  in  the 
following  sub  sections. 

3.1.1  High-Pressure  Air  Supply  and  Compressors 

The  storage  tank  volume  was  28-m3  with  a  maximum  allowable  design  pressure 
of  1.38  MPa.  The  Century  Electric  Co.  compressors  simultaneously  filled  up  the  storage 
tank  to  1.25  MPa.  The  compressors  were  run  by  two  20  hp  squirrel-case  induction  poly¬ 
phase  motor  (model  number:  SC-286U-FC3-3).  Both  compressors  at  full  running 
condition  usually  took  about  90  minutes  to  fill  up  the  tank.  An  Aerolab  heater  was  used 
to  dry  the  compressed  air  to  a  dew  point  of  -70  °C.  The  heater  had  two  tanks  filled  with 
silica-gel  desiccants.  The  compressed  air  passed  through  one  tank  while  the  desiccant  in 
the  other  tank  was  heated.  After  a  six-hour  cycle,  the  dryer  tanks  were  manually 
switched.  About  midway  through  the  project,  two  new  compressors  (Ingersoll-Rand 
model  SSR-HXP50SE)  were  purchased  at  no  additional  cost  to  the  ARO.  The 
compressors  rated  operating  pressure  was  1.38  MPa.  The  motor  used  for  each  of  the 
compressors  was  50  hp  (3  phase/  60  Hz).  An  automatic  air-cooled  Ingersoll-Rand  dryer, 
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model  TZ300HP-EMS-3V-LDP  (Desiccant-  activated  alumina,  dew  point  of  -73  °C,  was 
also  installed.  The  pump-up  time  with  the  new  compressors  was  about  30  minutes. 

3.1.2  Stagnation  Chamber 

The  stagnation  chamber  consisted  of  flow  straightening  tubes,  a  damping  screen 
and  a  Pitot  probe.  The  flow  straightening  tubes  were  10.16  cm  long  and  1.27  cm  in 
diameter.  The  steel  tubes  were  metal-jacketed  with  an  upstream-perforated  cone  and  a 
downstream  damping  screen.  A  schematic  drawing  of  the  flow  straightener  is  shown  in 
Fig.  4. 

3.1.3  Mach  3  Nozzle  &  Operating  Conditions 

The  Mach  3  nozzle  was  made  of  aluminium.  The  measured  Mach  number  was 
2.9710.04  (see  Appendix  A).  The  nozzle  block  was  25.4  cm  long,  with  a  1.65  x  7.62  cm2 
throat  area.  A  schematic  drawing  is  given  in  Fig.  5.  A  Pitot  tube  was  inserted  through  the 
nozzle  block  to  measure  the  upstream  total  pressure.  The  upstream  total  temperature  was 
monitored  with  an  Omega  brand  type  K  thermocouple  (the  uncertainty  was  ±  2K).  The 
tunnel  operating  pressure  was  regulated  by  the  compressed  air  (from  a  different  source) 
that  was  coupled  with  a  7.62  cm  (3  inch)  LESLIE  control  valve.  The  nominal  operating 
pressure  of  the  tunnel  was  0.64  ±  0.005  MPa  for  the  mean  flow  measurements  using  five- 
hole  probe.  For  the  PIV,  and  shadowgraph/schlieren  measurements  the  nominal  pressure 
was  0.61  ±  0.021  MPa.  The  tunnel  operating  conditions  are  summarized  in  Table  3. 

3.1.4  Test  Section 

The  dimensions  of  the  test  section  were  7.62  x  7.62  x  30.5  -cm3  [Fig.  6(a)].  The 
wind  tunnel  models  were  bolted  to  the  floor  of  the  test  section,  as  shown  in  Fig.  6(b).  The 
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sidewalls  contained  optical  glass  windows.  A  test  section  ceiling  made  of  solid  optical 
grade  Plexiglas  was  used  for  the  PIV,  PSP  and  surface  oil  flow  visualization.  A  slotted 
solid  aluminium  ceiling  was  used  for  the  mean  flow  measurement  using  the  five-hole 
Pitot/cone-static  probe.  The  ceiling  was  slotted  in  two  different  longitudinal  locations 
such  that  the  probe  tip  could  be  traversed  at  0.86  cm  and  at  6.35  cm  downstream  of  the 
injector  port.  The  five-hole  probe  was  attached  to  the  stepper  motor  slide  for  the  vertical 
traverse  and  the  motor  was  mounted  on  a  rail  for  horizontal  traverse.  The  axial  position 
uncertainty  was  ±0.8  mm,  and  the  spanwise  traverse  location  uncertainty  was  ±0.4  mm. 

3.1.5  Diffuser  and  Muffler 

The  variable  area  supersonic  diffuser  was  designed  to  operate  over  a  range  of 
Mach  numbers,  1 .5  -  4.0.  The  diffuser  design  was  a  part  of  an  undergrad  student  project, 
and  later  was  used  for  this  research  work.  The  inlet  and  the  exit  area  of  the  diffuser  were 
fixed  at  7.62  cm  *  7.62  cm  and  12.7  cm  x  7.62  cm,  respectively.  The  hinged  plate  was 
pinned  to  the  walls  and  could  be  adjusted  with  variable  Mach  number  supersonic  wind 
tunnel  test  sections.  Diffuser  dimensions  are  listed  and  shown  in  the  Table  4  and  Fig.  7.  A 
silencer  was  constructed  and  attached  to  the  wall. 


3.2  Missile  Injector  Models 

3.2. 1  Missile  Fuselage  and  Injectors 

The  wind  tunnel  models  were  half-body  3-D  blended-wall  missile  fuselages  (Fig. 
2b).  A  single  port  circular  injector  was  machined  into  each  model.  The  models  were 
19.05  cm  long,  1.78  cm  in  diameter  ( D ),  and  13°  half-cone  angles.  A  schematic  of  the 
models  can  be  seen  in  Fig.  8.  Four  models  with  different  injectors  were  used  for  the 
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current  research.  The  models  were  designated  M3A-HP  (Mach  3,  air  injectant,  high 
injection  pressure),  M3A-LP  (Mach  3,  air  injectant,  low  injection  pressure),  M3H  (Mach 
3,  helium  injectant),  and  MIA  (Mach  1,  air  injectant).  A  fifth  model  (Ml  A-L,  Mach  1,  air 
injectant,  and  Larger  injector  diameter  of  0.953  cm)  was  also  examined  in  few 
preliminary  experiments,  but  the  disturbances  were  too  large  for  the  present  experimental 
apparatus.  Figure  9  shows  a  photograph  of  all  the  models.  Model  M3H  was  operated  with 
compressed  Helium  gas  and  the  rest  were  with  compressed  air.  Detailed  injector 
dimensions  are  listed  in  Table  5  and  the  injector  operating  conditions  are  summarized  in 
Table  6.  The  injector  operating  conditions  as  well  as  the  free  stream  conditions  were  set 
to  achieve  a  range  of  low  MPRs,  0.005  to  0.09. 

3.2.2  Pebble  Bed  Heater 

A  pebble  bed  heater  was  designed  and  constructed  for  high  temperature  injection. 
The  heater  dimensions  are  shown  in  Fig.  10.  The  hot  jet  injection  surface  oil  flow 
visualization  was  performed  at  1190  °F  (916  K).  The  pebble  bed  heater  and  the  electric 
furnace  (Fisher  Scientific  Model-14,  Cat  No.  10-750-14  &  14A)  used  for  this  experiment 
are  shown  in  Fig.  11.  The  ceramic  pebbles  were  first  heated  in  the  electric  furnace  to 
2000  °F  and  then  poured  into  the  home-built  heater.  The  heater  was  then  placed  under  the 
tunnel  test  section  [Fig.  11(d)],  Due  to  safety  concerns,  this  hot  jet  interaction  flowfield 
visualization  was  performed  only  for  the  sonic  model.  The  temperature  of  the  air  passing 
through  the  hot  pebble  bed  heater  was  recorded  with  two  thermocouples.  One  of  the 
thermocouples  measured  the  hearth  temperature  of  the  pebble  bed  heater  and  the  other 
measured  the  exit  air  temperature  (i.e.;  the  jet  total  temperature).  As  suggested  in  Chapter 
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2,  helium  was  used  to  simulate  higher  temperatures  than  were  available  with  the  pebble 
bed.  The  available  temperature  operating  conditions  are  listed  in  Table  7. 


3.3  Measurement  Apparatus 

3.3.1  Surface  Oil  Flow  Visualization  (SOFV) 

Dow  Coming  200  fluid  (1000CST)  was  used  for  surface  oil  flow  visualization. 
The  oil,  mixed  with  orange  fluorescent  dye  (Day-Glo  color  corporation,  pigment  sample 
T-15  blaze  orange),  was  smeared  onto  the  test  surface  and  illuminated  by  an  ultraviolet 
lamp.  A  CCD  camera  (Kodak  Professional  DCS  460c,  resolution  3060x2036  with  a 
Nikon  N90-  50mm  lens)  mounted  on  top  of  the  test  section  was  used  to  collect  the  images 
of  the  flow  structure  during  the  tunnel  run  condition.  The  camera  was  attached  on  to  the 
tri-pod  and  images  were  taken  through  the  transparent  test-section  ceiling  made  of  optical 
grade  solid  Plexiglas.  The  camera  was  connected  (by  a  SCSI  port  and  a  cable)  to  an  NEC 
Intel  Pentium  II  PC  running  Microsoft  Windows  98.  An  Adaptec  AVA-1502  SCSI  Host 
Adapter  with  a  SCSI2  port  transferred  the  digital  image  recorded  by  the  CCD  imager  to 
the  computer.  The  image  was  then  acquired  and  displayed  through  the  Kodak  TWAIN 
driver  into  the  Adobe  Photoshop  4.0  application  software  for  Windows. 

3.3.2  Five-hole  Pitot  /Cone-Static  Pressure  Probe 

An  Aeroprobe  Inc.  Brand  L-shaped  five-hole  Pitot/cone-static  pressure  probe  with 
half  cone  angle  of  19°  was  used  to  measure  the  mean  Mach  number.  Figure  12  shows  the 
5-hole  probe.  A  stepper  motor  with  velmex  bi-slide  rail  was  used  to  traverse  the  probe  in 
the  vertical  y-direction.  The  traverse  distance  was  about  5.0  cm  from  the  test-section 
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floor.  A  Linear  Voltage  Displacement  Transducer  (LVDT),  Omega  model  LD100  with  a 
SP200A  controller  was  used  to  record  the  probe  height.  The  analog  output  of  the  LVDT 
transducer  was  set  to  0.0  -  5.0  volts  for  a  range  of  movement  of  0.0  -  5.0  cm.  The  x  and  z 
locations  were  manually  set. 

3.3.3  Two-Color  Digital  Particle  Image  Velocimetry  (PIV) 

The  digital  two-color  PIV  system  consisted  of  two  Continuum  brand  Surelite 
Nd:YAG  lasers  that  were  powered  by  two  Continuum  Model  SLI-10  lasers.  The  surelite 
lasers  each  produce  coherent  green  (532  nm,  225  mJ)  laser  beams.  One  of  the  green  laser 
beams  was  routed  through  a  red  dye  laser  (607  nm).  A  12-volt  DC  motor  driven  gear 
pump  was  used  to  circulate  the  dye  (sulfarhodamin  640  dissolved  in  methanol).  The  dye 
laser  was  about  20-30%  efficient.  The  red  and  green  beams  were  then  co-located  in  space 
and  converted  into  a  single  sheet  by  a  series  of  optical  components.  The  time  delay 
between  the  red  and  green  beams  was  controlled  with  an  ISSI,  Inc.  custom  control  circuit. 
Model  DG535  four-channel  digital  delay/pulse  generator  (included  with  the  complete 
PIV  system)  was  used  to  set  a  time  delay  of  400  ns.  The  optics  that  directed  the  beams 
into  the  test  section  were  mounted  on  to  the  nearby  wall  to  avoid  vibration  during  tunnel 
runs.  A  Kodak  model  460c  CCD  camera  with  a  Nikon  N90  (150  mm)  lens  was  placed 
apart  from  the  wind  tunnel  table  (base)  to  prevent  vibration.  The  CCD  camera  captured 
the  image  through  optical  grade  glass  windows  in  the  test  section  sidewalls.  Because  of 
safety  concerns  with  the  high-pressure  seeding;  only  the  sonic  missile  models  were  used 
for  the  PIV  measurement. 

Atomized  olive  oil  was  used  as  the  seed  material.  A  TSI  brand  model  9306  six-jet 
atomizer  was  used.  The  atomizer  was  placed  in  a  compressed  air  chamber  [Fig.  13(a)]  to 
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ensure  proper  seeding  capability  at  higher  than  recommended  pressures.  The  flow 
tracking  capability,  3-dB  point  frequency  response  of  60-200  kHz  is  described  in 
Bowersox  (1996).  In-house  software  (b-PIVbeta-1.08)  was  used  to  perform  the  post 
processing  of  the  data  and  to  calculate  the  instantaneous  velocity  distribution. 

Due  to  camera  framing  rate  limitations  (1-2  images  per  tunnel  run),  this  system 
was  not  suitable  for  turbulence  measurements.  Hence,  PIV  measurements  were  acquired 
with  a  newer  alternate  PIV  system,  as  described  in  the  next  section. 

3.3.4  Lavision  Flowmaster3  PIV  System  for  2-D  Turbulence  Measurement 

The  LaVision  PIV  system  uses  a  Flowmaster3  CCD  camera  with  a  New  Wave 
Research  Nd-Yag  laser.  The  1280  x  1024  pixel,  12-bit  CCD  camera  provided  a 
maximum  of  8  single  frame  images  per  second.  Detailed  camera  specifications  are 
documented  in  Table  8.  The  532  nm  class  IV  laser  (detailed  specifications  are  listed  in 
Table  9)  and  the  camera  interface  were  all  linked  together  with  an  externally  triggered 
Programmable  Timing  Unit  (PTU)  and  a  PC.  Figure  13(b)  shows  the  components  of 
LaVision  PIV  system.  The  laser  power  was  set  at  60%  with  the  time  delay  between  the 
two  exposures  of  600ns.  Double  frame/  double  exposure  cross  correlation  mode  was  used 
to  capture  images.  The  focal  point  of  the  laser  sheet  was  placed  on  top  of  the 
measurement  plane  surface  that  provided  about  a  maximum  of  1.0  mm  sheet  thickness 
throughout.  The  PIV  optics  and  the  camera  were  mounted  onto  the  lab  wall  and  floor, 
respectively,  so  that  the  minimal  tunnel  vibration  would  not  affect  the  measurements. 
Olive  oil  was  seeded  with  a  TSI  Model  9306  six-jet  atomizer.  The  atomizer  (seeder)  was 
placed  in  a  tank  as  described  in  Section  3.3.3.  When  the  smoke-full  tank  reached  the 
required  pressure,  the  tank  outlet  valve  was  released,  and  thus  allowed  the  smoke  to  pass 
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through  the  injector  port.  Davis6.04  software  and  the  vector  analysis  3.0  software  were 
used  to  post  process  the  images. 

3.3.5  Pressure  Sensitive  Paint  (PSP) 

ISSI,  Inc.  brand  Uni-FIB  Pressure  Sensitive  Paint  (PSP)  was  used.  The  paint 
contained  two  parts:  an  oxygen  permeable  flouroacrylic  co-polymer  (FIB)  binder  and  a 
fluorinated  Platinum  porphyrin  meso-Tetra-pentaflourophenyl  porphrine  [Pt(TfPP)] 
oxygen  sensitive  probe  molecule.  The  PSP  was  applied  using  an  airbrush  sprayer.  The 
temperature  sensitivity  was  0.6%/°C.  The  paint  was  illuminated  with  two  blue  (464  nm) 
LED  light  sources  (0.25  W  each).  The  luminescence  output  (650  nm)  was  acquired  with  a 
Pixelvision  SpectraVideo  16-bit  CCD  camera  (model  SV512V1A/PFT-95).  The  paint 
calibration  stemmed  from  the  Stem-Volmer  relationship,  where  the  pressure  was 
expressed  as  a  function  of  Ire/I.  As  suggested  by  the  manufacturer,  the  calibration  data 
were  fit  with  a  second  order  polynomial.  The  uncertainty  for  the  PSP,  as  provided  by  the 
manufacturer  was  ±2.0%  of  the  measured  value.  The  overall  uncertainty  is  described  in 
Appendix  B. 

3.3.6  Schlieren/  Shadowgraph  Optics 

The  schlieren/shadowgraph  system  consisted  of  two  concave  mirrors  (20.32cm 
Dia,  99.06cm  FL),  one  light  source,  a  screen,  a  CCD  camera,  a  color  knife-edge,  and  a 
PC.  Shadowgraph  setup  was  relatively  simple.  A  single  concave  collimating  mirror  was 
used  to  create  a  collimated  light  beam  that  passed  through  the  test  section  sidewall 
windows  and  onto  a  screen.  The  Kodak  460dcs  camera  system  (described  above)  was 
used  to  capture  the  image  on  the  screen.  For  the  schlieren  images,  the  collimated  beam 
that  passes  through  the  test  section  was  projected  onto  a  second  mirror.  The  color  knife- 
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edge  was  located  at  the  focal  point  of  the  reflected  beam.  The  diverging  beam  was 
viewed  on  the  screen,  and  the  CCD  camera  captured  the  color  image. 

3.3. 7  Pressure  Transducers 

A  number  of  pressure  transducers  were  employed  for  the  pressure  measurements 
at  different  locations.  For  the  upstream  total  pressure  and  the  test-section  static  pressure 
measurements  Omega  brand  model  PX302  (range:  0-100  psi)  and  PX  202  (range:  0-15 
psi)  were  used.  DP25  signal  conditioners  were  used  for  both  of  these  transducers.  For  the 
mean  flow  measurement  using  the  five-hole  Pitot/cone-static  pressure  probe  four  0-1 5psi- 
range  Omega  model  PX139-015A4V  and  one  0-100psi-range  omega  model  PX139- 
100A4V  transducers  were  used.  Transducer  specifications  are  given  in  Table  10. 
Calibration  of  each  transducer  and  the  detailed  uncertainty  analysis  are  described  in 
Appendix  B. 

3.3.8  Traverse  System 

The  traverse  system  consisted  of  a  Linear  Voltage  Displacement  Transducer 
(LVDT)  with  a  signal  conditioner,  a  stepper  motor,  a  motor  controller,  a  velmex  bi-slide, 
mounting  to  hold  the  probe,  and  the  horizontal  rail  for  the  manual  horizontal  traverse. 
Omega  Model  LD 100- 150  LVDT  with  nominal  linear  range  152.4mm  was  used  to  record 
the  vertical  position  of  the  probe.  The  LVDT  consisted  of  a  precision  variable  inductor 
with  frequency  response  up  to  15khz,  and  an  aluminum  core  that  moved  in  a  polyimide 
tube.  The  core  movement  transmitted  equivalent  amount  of  inductance  to  the  signal 
processor.  Omega  Model  SP200A.  The  raw  voltage  produced  from  the  signal  processor 
was  linearly  proportional  to  the  inductance.  The  uncertainty  of  the  LVDT  was  0.07  mm. 
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3.3.9  High  Temperature  Thermocouple 


For  high  temperature  surface  oil  flow  visualization,  the  pebble  bed  heater  was 
equipped  with  two  Omega  model  XTA-W5R26-U-125-30-H-HX-6  thermocouples.  The 
maximum  operating  temperature  of  each  thermocouple  was  2300°C.  Thermocouple 
uncertainties,  as  specified  by  the  manufacturer,  were  0.4  °C. 

3.3.10  Data  Acquisition  System 

The  data  acquisition  system  contained  an  Analog  to  Digital  (A/D)  conversion 
card,  an  A/D  board  (16-channel  TSI  model  ADCWIN-16),  a  Gateway2000  pentium  PC, 
and  the  TestPoint  software  provided  by  the  Keithley  Instruments.  A  program  in  the 
TestPoint  software  environment  was  written  to  record  the  pressure  and  temperature  data. 
For  SOFV  and  PIV,  the  static  and  the  upstream  total  pressure  measurements  were  taken 
at  a  sampling  rate  of  1000  Hz  for  a  total  of  12000  points  (tunnel  runtime  of  12  seconds). 
For  the  five-hole  Pitot/cone-static  mean  flow  measurement,  a  sampling  rate  of  1000  Hz 
was  set  for  a  total  of  12000  points  during  the  vertical  traverse  from  on  top  of  the  missile 
model,  and  a  total  of  15000  points  were  set  while  traversing  from  on  top  of  the  test 
section  floor. 
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4.  EXPERIMENTAL  METHODOLOGIES 


4.1  Five-Hole  Probe 


The  upstream  total  pressure  Pt00,  four  cone  surface  pressures  (Pci,  PC2,  PC3  and  PC4), 
and  Pitot  pressure  (Pt2)  were  recorded.  Two  axial  locations  were  selected  for  the  mean 
flow  measurement;  one  at  15.34  mm  downstream  of  the  injector  center  and  the  other  was 
at  46.96  mm  downstream  in  the  wake  region.  At  each  axial  location,  the  probe  was 
traversed  spanwise  at  15  locations.  The  spanwise  locations  at  the  center  of  the  test  section 
were  0.318  cm  apart.  A  FORTRAN  program  was  written  to  post-process  these  raw 
voltage  data.  The  program  first  averaged  and  calibrated  the  LVDT  data.  The  program 
produced  180  data  points  across  the  height  of  about  4.7  cm  from  the  test-section  floor, 
where  each  point  consisted  of  an  average  of  80  data  samples. 

Absolute  pressure  data  were  obtained  using  the  transducer  calibration  equations 
(listed  in  Appendix  B).  For  a  range  of  different  supersonic  Mach  numbers  spanning  that 
initially  expected  in  the  present  experiments,  the  ratio  of  cone  surface  pressure  (Pc)  and 
Pitot  pressure  (Pt2)  was  obtained  from  the  numerical  solution  of  Taylor-Maccoll  conical 
flow  equation.  The  result  is  given  in  Table  11.  From  this  theoretical  result  a  least  squares 
curve  fit  solution  for  the  Mach  number  as  a  function  of  Pc/Pt2  was  obtained  as 
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The  Taylor-Maccoll  conical  flow  solution  was  valid  for  attached  shock,  but  the 
five-hole  probe  tip  was  blunt,  so  a  detached  shock  wave  was  produced.  Thus,  the 
theoretical  solution  was  not  directly  applicable.  To  minimize  errors,  experimental  results 
were  combined  with  the  theoretical  results  and  plotted  in  Fig.  14  for  a  range  of  available 
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Mach  numbers  1.72,  3,  4,  5.0  and  5.5.  Based  on  these  data,  a  refined  calibration  equation 
was  found  as 

M  =  (-0.0606  M,  +  1 .2299  )M,  (4.2) 

A  comparative  study  on  the  results  from  the  conical  theory  solution,  normal  shock 
solution  and  the  recalibrated  solution  are  shown  in  Table  12. 

A  FORTRAN  program  later  computed  other  local  mean  flow  quantities  using  the 
usual  compressible  flow  relations  [Anderson  (1990)].  The  Mach  number  contour,  the 
upstream  total  pressure  ratio  (Pti/Ptoo)  contour,  and  the  downstream  total  pressure  ratio 
(Pt2/Pt<»)  contour  were  plotted  with  Tecplot  version  8.0.  The  upstream  total  pressure  Pti 
was  found  from  the  fundamental  shock-wave  relationships.  The  uncertainty  associated 
with  the  Mach  number  measurement  by  the  five-hole  probe  was  estimated  as  ±0.13.  The 
uncertainty  calculations  are  described  in  Appendix  B. 

4.2  Pressure  Sensitive  Paint  (PSP) 

Pyrene  based  (PtOEP)  pressure  sensitive  paint  was  used  to  map  the  surface 
pressure  distribution  over  the  missile  model.  The  Pixelvision  SpectraVideo  CCD  camera 
was  used  to  capture  the  background,  wind-off  and  wind-on  images.  For  the  12-second 
tunnel  run  time,  five  images  were  taken.  ISSI  Inc.  PSP_Analysis_2  software  was  used  to 
average  the  images.  Averaged  background  image  was  then  subtracted  from  both  the 
averaged  wind-off  and  averaged  wind-on  images.  The  ratio  of  the  subtracted  wind-off 
and  the  subtracted  wind-on  images  was  then  converted  into  the  surface  pressure.  The 
PSP_Analysis_2  was  also  equipped  with  the  calibration  file,  so  that  the  software  directly 
converted  the  ratio  file  into  the  surface  pressure  (psiA)  using  the  estimated  surface 


24 


temperature  of  about  10  °C  (±5.0°).  The  overall  uncertainty  was  3.6%  of  the  measured 
value  (see  Appendix  B). 


4.3  La  Vision  PIV 

The  LaVision  Inc  flowmaster3  PIV  system  was  used  to  obtain  axial  and 
transverse  velocity  turbulence  information.  The  measurements  included  longitudinal 
contours  of  the  mean  velocity  field,  the  mean  vorticity  field,  the  mean  shear  streain  rate 
field,  velocity  component  turbulence  intensities,  turbulent  stresses,  Turbulent  Kinetic 
Energy  (TKE)  and  the  mean  dissipation.  For  each  model  tested,  x-y  plane  contours  were 
acquired  at  z  =  0,  1.59mm,  and  3.18mm.  Davis  6.04  software  (integrated  with  the  PIV 
system)  was  used  to  acquire  and  process  the  Mie-scattering  images.  Due  to  the  small 
lateral  extent  of  the  plumes,  no  particles  were  present  in  the  plane  located  at  z  =  4.76  mm. 
The  sonic  model  (model  MIA  at  atmospheric  pressure)  and  the  Mach  3  model  (model 
M3A-LP  at  pressures,  377  kPa  and  694.3  kPa)  were  investigated  using  this  PIV  system. 
At  each  measurement  location  about  30  tunnel-runs  were  performed,  where  each  run 
provided  40  instantaneous  velocity  field  measurements.  Due  to  continuous  seeding,  the 
smoke  particles  accumulated  and  liquefied  inside  the  tubing.  Images  that  were  splashed 
with  seeder  oil  were  unusable.  During  image  processing  the  contaminated  images  were 
deleted. 

Davis6.04  software  was  used  to  obtain  the  instantaneous  velocity  field 
information  at  each  run.  Cross  correlation  mode  was  used  to  process  the  PIV  images. 
Standard  11*12  [via  FFT,  no  zero  padding]  correlation  function  with  32  x  32  interrogation 
window  size  and  50%  overlap  were  employed  to  obtain  the  vector  field.  The  detailed 
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parameters  that  were  set  to  process  the  raw  images  and  to  post  process  the  vectors  are 
presented  in  Table  13.  The  vector  post  processing  was  done  in  two  steps.  The  filters  were 
first  applied  to  the  instantaneous  vector  field,  and  then  the  vector  length  cut-off  was  used 
to  obtain  the  correct  instantaneous  vectors.  The  valid  vectors  were  then  saved  in  ASCII 
data  file  format.  Each  data  file  consisted  of  5120  (80  x  64)  data  points,  corresponding  to 
an  image  area  of  about  67  mm  x  53  mm.  The  data  files  were  then  used  to  obtain  the 
turbulence  information  using  an  in-house  FORTRAN  program  that  was  written  for  this 
purpose. 

The  FORTRAN  program  first  computed  the  average  velocity  at  each  grid  point. 
The  average  was  computed  if  the  total  number  of  valid  vectors  at  a  point  exceeded  the 
number  of  at  least  70%  of  the  total  number  of  possible  vectors  at  that  point.  To  estimate 
the  gradients  in  the  transport  equations  a  first  order  central  difference  approximation  was 
used  within  the  plume  area  and  the  forward  and  backward  difference  approximations 
were  used  at  the  plume  upper  and  lower  edges,  respectively.  The  turbulence  quantities 
were  estimated  using  the  following  equations: 
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The  freestream  velocity  and  missile  diameter  (D 
nondimensionalize  the  turbulence  results. 


5.  FLOWFIELD  RESULTS 


The  injector  operating  conditions  are  summarized  in  Tables  1,  6  and  7.  The  tunnel 
operating  conditions  are  summarized  in  Tables  2  and  3.  As  indicated  in  Table  6,  a  wide 
range  of  injector  conditions  was  examined.  Of  particular  importance  to  the  subsequent 
results  is  the  effective  backpressure  ratio  [Schetz  and  Billig  (1967)]  that  varied  over  a 
range  of  0.19-6.15,  which  indicated  that  underexpanded  (pefpx  >  1),  perfectly  expanded 
(Pei/pcc  =  1)  and  overexpanded  (pei/poo  <  1)  injection  conditions  were  present.  The 
majority  of  the  low-blowing  ratio  cases  studied  here  were  overexpanded. 

5.1  Flow  Visualization 

5.1.1  Surface  Oil  Flow  Visualization  (SOFV) 

The  flow  surface  interaction  topologies  were  discerned  from  the  Surface  Oil  Flow 
Visualization  (SOFV)  experiments.  The  salient  flow  features  include  (1)  the  missile  body 
shock  at  the  tip  of  the  missile  fuselage,  (2)  the  boundary  layer  separation,  (3)  the  jet 
interaction  region,  (4)  the  interaction  shock,  (5)  the  horseshoe  vortex,  (6)  the  re¬ 
compression  shock  and  (7)  the  wake  region.  These  features  are  annotated  in  Fig  15  for 
the  sonic  injection  case  at  an  MPR  of  0.006.  For  this  case,  the  jet  exit  is  nearly  perfectly 
expanded  (effective  back  pressure  ratios  listed  in  Table  6  was  very  close  to  unity).  In  Fig. 
15,  the  freestream  flow  is  directed  from  left  to  right  and  the  jet  flow  is  directed  outward 
from  the  paper.  Flowfield  symmetry  is  clearly  depicted  from  the  image.  The  missile-body 
shock  angle  of  24°  was  estimated  from  the  image,  while  the  shock  angle  calculated  from 
conical  theory  was  23.8°.  When  the  crossflow,  downstream  of  the  missile  shock, 
encountered  the  expansion  edge  of  the  missile  fuselage,  it  deflected  and  marched  parallel 
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to  the  cylindrical  missile  fuselage.  The  jet  plume  acted  as  an  obstacle  to  the  crossflow; 
this  resulted  in  the  boundary  layer  separation  upstream  of  the  jet  and  created  the 
interaction  shock.  The  round  thick  line  of  oil  accumulation  represents  the  separation  line 
due  to  the  interaction  shock.  The  lambda  shock  in  between  the  boundary  layer  separation 
point  and  the  interaction  shock  cannot  be  identified  in  SOFV  images.  However,  in 
shadowgraph  and  schlieren  images  (discussed  later),  the  lambda  shock  is  clearly 
distinguished.  The  next  thin  oil  line  that  wrapped  around  the  jet  was  interpreted  as  the 
horseshoe  vortex  [Everett  et  al  (1998)].  Two  other  distinct  oil  lines  were  observed 
downstream  of  the  jet.  These  lines  were  created  due  to  the  wake  vortices  [Fric  and 
Roshko  (1994)].  The  crossflow  stream  that  wrapped  around  the  jet  merged  to  create  the 
compression  shocks  that  appeared  downstream  of  the  jet. 

A  comparison  of  the  SOFVs  for  model  MIA  is  presented  in  Figs.  16(a  -  c).  The  flow 
features  for  this  model  at  MPR  0.006  have  already  been  discussed  above;  the  jet  total 
pressure  used  for  this  case  was  101.35  kPa.  Figure  16(b)  shows  the  same  model  at  an 
MPR  of  0.005.  This  reduced  MPR  was  achieved  by  operating  the  injector  at  a  lower  jet 
total  pressure  of  87.56  kPa.  Although  the  jet  flow  in  Fig.  16(b)  was  overexpanded  (see 
Table  6),  the  surface  flow  features  were  very  similar  to  those  for  an  MPR  of  0.006  [Fig. 
16(a)].  Figure  16(c)  shows  the  high  temperature  injection  case.  Again,  the  flow  features 
are  very  similar  to  those  in  Fig.  16(a).  Champigny  and  Lacau  (1994)  reported  that  higher 
jet  temperatures  shift  the  separation  line  further  upstream,  however  in  the  present  study,  a 
300%  increase  in  the  injection  static  temperature  resulted  in  only  a  4.0%  increase  in  the 
boundary  layer  separation  distance.  The  most  likely  cause  for  this  discrepancy  is  that  the 
low-blowing  ratio  data  in  Fig.  16  are  either  perfectly  expanded  or  slightly  overexpanded, 
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and  most  of  the  cases  presented  in  the  literature  are  for  high  blowing  ratios,  which  are 
highly  underexpanded.  Underexpanded  jets  create  larger  upstream  flow  disturbances  than 
overexpanded  jets,  thus  the  temperature  effects  are  likely  to  be  more  significant.  The 
separation  distance  is  discussed  in  more  detail  later  in  this  section. 

The  SOFV  for  model  M3A-LP  is  presented  in  Fig.  17.  Figure  17(a)  corresponds  to 
an  MPR  of  0.03  and  Fig.  17  (b)  to  0.06.  For  these  two  cases,  the  injector  was  significantly 
overexpanded,  (see  Table  6).  The  flow  features  are  analogous  to  those  seen  for  sonic 
injection  case.  The  primary  differences  were  observed  are  in  the  wake  region. 
Specifically,  the  oil  lines  in  Fig.  17(a)  immediately  downstream  of  the  jet  moved  further 
apart  while  progressing  downstream  until  about  2.5 De  from  the  injection  port.  Then,  the 
two  lines  merged  (~  4.5 De  downstream  of  the  exit)  creating  a  diamond  shaped  structure. 

The  SOFV  for  model  M3A-HP  is  given  in  Figure  1 8.  High  jet-to-freestream  pressure 
ratios  of  9.4  and  18.8  were  used  for  this  Mach  3  injector,  which  provided  momentum 
ratios  (MPR)  of  0.05  and  0.09,  respectively.  Here,  the  jet  flows  were  underexpanded. 
Aside  from  separation  distance,  the  differences  between  the  low  and  high  MPR  images 
were  only  discernible  in  the  wake  region.  For  the  higher  MPR,  the  jet  penetration  was 
higher,  and  the  wake  region  developed  further  downstream.  This  is  evident  from  the 
wake  lines  that  moved  downstream  at  small  angle  until  about  10De  downstream  of  the 
injector  port  [Fig  18(b)],  and  then  the  angle  expanded  outward  indicating  a  developed 
wake  region.  However,  for  the  low  MPR  case  [Fig  18(a)],  the  small  angle  prevailed  until 
about  lDe  downstream. 

Like  for  the  M3A-LP  model,  two  momentum  ratios,  0.03  and  0.05,  were  examined 
for  the  Mach  3.0  helium  injection  with  model  M3H.  The  results  are  presented  in  Figs 
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19(a  and  b),  respectively.  The  SOFVs  for  model  M3H  in  Fig.  19,  strongly  resembled 
those  for  model  M3A-LP,  as  shown  in  Fig  17,  where  the  injector  operating  conditions 
were  matched.  For  this  case,  the  use  of  helium  resulted  in  an  equivalent  air  injection 
temperature,  based  on  velocity  ratio,  of  1250  K. 

To  simulate  an  even  higher  temperature  air  injection,  the  MIA  injector  was  operated 
at  MPR  0.015  with  helium.  Air  was  also  used  for  comparative  purposes.  The  use  of 
helium  simulated  air  injection  at  a  temperature  of  1900  K.  Again,  the  injector  operating 
conditions  are  summarized  in  Tables  6  and  7.  Figures  19(c  -  d)  present  the  air  and  helium 
SOFV,  respectively.  The  flow  structures  looked  similar.  The  impact  of  temperature  on 
separation  distance  is  discussed  in  the  next  paragraph. 

The  separate  distances,  as  inferred  from  the  SOFVs  in  Figs.  16-19,  are  summarized 
in  Table  14.  The  uncertainties  are  described  in  Appendix  B.  Figure  20  presents  the 
variation  of  the  separation  distance,  normalized  by  injector  diameter,  with  respect  to  the 
MPR,  for  all  of  the  test  cases.  Overall,  the  expected  trend  of  increasing  of  the  separation 
distance  with  increasing  MPR  was  observed.  However,  MPR  alone  is  not  a  sufficient 
correlation  parameter,  where  systematic  increases  in  separation  distance  with  increasing 
injection  pressure  and  temperature  were  discerned.  For  example,  the  M3A-HP  separation 
distances  were  significantly  above  the  M3A-LP  trend.  An  explanation  for  this  trend  can 
be  found  by  referring  back  to  Table  6.  Specifically,  the  M3A-HP  cases  were 
underexpanded,  and  hence  the  plume  expanded  upon  entering  crossflow.  Thus,  the  size  of 
the  obstruction  increased,  as  compared  to  the  M3A-LP  case,  where  the  flow  was 
overexpanded.  Figure  21  shows  the  variation  in  the  separation  distance  (normalized  by 
De )  with  the  jet-to-freestream  pressure  ratio  for  the  M3A  models.  Figure  22  shows  the 
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effect  of  injection  temperature  on  separation  distance.  Over  the  range  of  conditions  tested 
here,  the  separation  distance,  for  a  given  MPR,  increased  at  a  rate  of  approximately  0.013 
with  the  jet-to-freestream  temperature  ratio.  It  is  anticipated  that  if  the  high  temperature 
jets  were  operated  in  an  over  expanded  condition,  the  influence  of  temperature  would 
have  been  more  substantial. 

5.1.2  Mie- Scattering  Visualization 

Mie  scattering  images  were  obtained  for  models  MIA  (MPR  =  0.006)  and  M3A-LP 
(MPR  =  0.03  and  0.06).  Mie-scattering  images  were  taken  using  both  the  ISSI  and 
LaVision  brand  PIV  systems.  Due  to  the  limitations  discussed  earlier,  the  former  system 
was  used  only  for  the  sonic  model  (MIA)  for  the  flow  visualization  purposes. 

Depicted  in  Fig.  23  is  the  Mie-scattering  image  of  model  Ml  A  for  the  MPR  of  0.006 
obtained  with  the  digital  two-color  PIV  system.  Figures  24  (a  -  c)  show  the  Mie- 
scattering  images  for  the  same  model  with  the  LaVision  PIV  system  at  the  three  spanwise 
locations.  The  appearance  of  the  plume  was  nearly  identical  as  acquired  by  the  two 
systems  [comparing  Fig.  23  to  Fig.  24(a)].  Hence,  the  subsequent  discussion  will  be  in 
reference  to  the  images  obtained  with  the  LaVision  system.  At  about  6.5De  downstream 
(visual  inspection)  of  the  injector  port,  the  jet-plume  surface  interaction  is  clearly 
observed.  Large-scale  turbulence  structures  are  clearly  visible  in  these  images.  Figure  25 
(a  -  c)  presents  the  time-averaged  Mie-scattering  images  for  the  sonic  model  at  the  three 
measurement  planes.  These  data  show  the  mean  penetration  and  where  the  plume 
interacted  with  the  missile  fuselage. 

The  instantaneous  and  average  images  for  the  model  M3A-LP  are  presented  in  Figs 
26  -  27  (a  -  c)  for  an  MPR  of  0.06  at  the  three  spanwise  planes.  The  plume,  upon  exiting 
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the  injector  port,  was  turned  downstream  by  an  oblique  shock  cone.  The  terminating 
oblique  shock  system  was  also  observed  in  the  shadowgraph  and  schlieren  images  (next 
section).  Figures  28  -  29(a  -  c)  present  the  instantaneous  and  average  Mie-scattering 
images  for  the  model  M3A-LP  at  lower  MPR  of  0.03.  Due  to  lower  MPR,  the  plume 
penetration  was  lower  than  that  for  the  higher  0.06  MPR,  which  resulted  earlier  jet- 
surface  interaction.  For  this  case,  the  terminating  oblique  shock  was  almost  normal  to  the 
jet  exit  port;  this  was  also  observed  in  shadowgraph  images  discussed  in  the  next  section. 
The  importance  here  is  that  the  overexpanded  flow  structure  turned  the  flow  downstream, 
and  thus  reduced  the  angle  and  size  of  the  jet  plume  induced  flow  blockage.  This  explains 
the  above  observed  difference  in  the  shock  separation  distance  when  comparing  the 
underexpanded  and  overexpanded  cases. 

5.1.3  Schlieren/  Shadowgraph  Images 

Shadowgraph  and  schlieren  images  were  acquired  for  qualitative  assessment  of 
the  shock  structures  for  the  different  injectors.  The  salient  features  in  the  images  include 
(1)  the  structure  of  the  missile  body  shock,  (2)  the  interaction  shock,  (3)  the  lambda 
shock,  (4)  the  recompression  shock,  and  (5)  the  exit  port  shock  structure  (barrel  shock  for 
underexpanded  injection  and  the  shock  cone  for  the  overexpanded  condition).  Also 
noticeable  were  weak  waves  (close  to  the  freestream  Mach  angle)  that  were  generated 
due  to  the  surface  imperfections  in  the  nozzle  and  test  section.  The  shadowgraph  and  the 
schlieren  images  are  presented  in  Fig  30(a  -  j)  and  Fig.  31  (a  -  f),  respectively.  The  flow 
features  (annotated  in  the  figures)  are  identical  in  both  the  shadowgraph  and  schlieren 
images.  Both  air  and  helium  injection  models  (MIA,  M1A-L,  M3A-HP,  M3A-LP,  and 
M3H)  were  examined  using  the  shadowgraph  technique.  Because  of  the  strong  similarity 


33 


between  the  shadowgraph  and  schlieren  images,  schlieren  photographs  were  only 
acquired  for  the  MIA,  M3A-LP,  and  M3A-HP  models.  In  Figs.  30  and  31,  the  freestream 
flow  was  from  right  to  left,  and  the  jet  flow  direction  was  vertically  upward  in  the 
positive  y-direction.  The  wave  angles  at  the  nozzle  exit  in  the  images  (19.5-20.0  degrees) 
agreed  well  with  the  Mach  wave  angle  of  19.7  degrees.  The  left  images  in  each  row 
correspond  to  the  lower  MPR  case.  With  the  exception  of  model  M3A-HP,  jet  plume 
barrel-shocks  were  not  visible  in  the  images.  Instead,  the  plume  was  terminated  by  the 
oblique  conical  shock  system.  A  barrel-shock  and  Mach  disk  are  faintly  visible  in  the 
high-pressure  model  M3A-HP  schlieren  [Fig.  31(f)]]. 

5.2  Mean  Flow  Measurements 

5.2.1  Five-hole  Probe 

The  five-hole  probe  was  used  to  measure  the  Mach  number  contours,  and  the 
pressure  ratio  contours  (Pa/Pti,  Pe/Pt®,  Pc/Pt®,  and  Pti/Pt®)  in  two  transverse  locations, 
15.34  mm  and  46.96  mm  downstream  of  the  injector  port.  The  coordinate  axes  and  key 
flow  features  are  labeled  in  the  Mach  number  contours  of  model  MIA  [Fig  32(a)],  The 
contour  plots,  for  both  the  upstream  and  downstream  measurement  locations,  are  drawn 
in  the  same  coordinate  axes  to  better  display  the  flowfield  evolution.  The  freestream  flow 
in  all  of  these  contour  plots  is  from  left  to  right;  the  jet  flow  is  in  the  positive  y-direction 
with  the  injection  port  centered  at  (0,0,0);  the  co-ordinate  axes  are  labeled  in  inches.  The 
M3H  model  was  not  examined  using  five-hole  probe  because  the  presence  of  helium 
contaminates  the  probe  Mach  number  calculation. 
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Presented  in  Figures  32  are  the  Mach  number  contour  plots  for  models  MIA 
(MPR  =  0.006,  0.005),  M3A-LP  (MPR  =  0.03,  0.06),  and  M3A-HP  (0.05,  0.06);  these 
cases  correspond  to  nearly  perfectly  expanded,  underexpanded  and  overpanded, 
respectively.  Figures  33  -  36  present  pressure  ratio  contour  plots  ( PJPti ,  Pti/Pua  Pa/Ptm 
and  P c/P ioa  respectively)  of  models  MIA,  M3A-LP,  and  M3A-HP.  At  the  upstream 
location,  the  flow  features  are  very  distinct  and  uniform.  The  general  flow  features  are 
annotated  in  Fig.  32  (a),  which  presents  Mach  number  contour  of  model  MIA  at  MPR 
0.006.  The  wake  region,  jet  plume  area,  and  the  interaction  shock  are  labeled  in  the 
figure.  The  flow  features  are  identical  in  the  pressure  ratio  plots.  The  plume  and  the 
interaction  shock  locations  were  verified  with  the  Mie-scattering  and  shadowgraph 
images.  The  horseshoe  vortices  at  both  sides  of  the  jet  plume  near  the  missile  surface 
were  also  noticeable.  A  small  region  in  the  wake  at  the  upstream  station  was  manually  set 
to  zero  because  of  flow  reversal  that  caused  lower  Pitot  pressures  (P^)  than  the  cone 
static  pressures  (pci,  pC2,  pC3,  or  pC4),  which  in  turn  resulted  erroneous  Mach  number  and 
pressure  ratio  estimations.  Flowever,  Figs.  35-36  present  the  five-hole  probe 
measurements  (Pitot  and  cone-static  pressures,  Pa/P  too  and  Pc/PvP)  without  any 
modification.  The  pressure  ratios  as  normalized  by  the  freestream  total  pressure, 
eliminate  the  errors  associated  with  the  stagnation  pressure  control  deviations. 

The  contour  plots  for  model  MIA  at  MPR  0.005  and  model  M3A-LP  at  MPR  of 
0.03,  demonstrated  similar  flow  features  to  those  as  described  above.  Presented  in  Fig  32 
(d)  is  model  M3A-LP  at  MPR  0.06,  where  the  flow  features  were  similar  to  those  as  seen 
in  the  sonic  injection  case,  except  the  interaction  shock  was  at  an  elevated  location,  which 
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indicated  higher  penetration  for  the  Mach  3  injection  case.  At  the  downstream  location, 
the  wake  merged  with  the  plume  to  create  a  large  circular  region  of  low  Mach  number. 

Mach  contours  for  model  M3A-HP  can  be  seen  in  Figs  32  (e  -  f).  Due  to  the  high- 
pressure  injection,  the  penetration  height  increased  and  the  interaction  shocks  were 
observed  only  at  the  upstream  location.  In  the  downstream  contour,  the  plume  shape  was 
not  clearly  discemable,  and  the  wake  region  appeared  to  be  the  most  prominent  structure. 
This  result  was  evident  for  the  higher  penetration  of  the  plume.  The  low  MPR  model 
M3A-LP  mean  flow  was  more  complicated  in  the  downstream  location.  This  was  because 
the  jet  remained  confined  close  to  the  surface,  which  caused  irregular  and  complicated 
surface  interactions. 

5.2.2  Pressure  Sensitive  Paint  (PSP) 

The  pressure  sensitive  paint  contours  for  the  same  cases  discussed  above,  with  the 
addition  of  the  M3H  case,  are  shown  in  Figs.  37-40.  The  freestream  flows  from  top  to 
bottom  in  all  the  images.  A  common  scale  with  a  range  from  0  to  6  psiA  was  employed 
for  all  models,  so  that  direct  comparisons  could  be  made.  The  entire  painted  image  area 
was  35.56  mm  x  83.35  mm.  This  region  covered  about  29.7  mm  upstream  and  about  53.7 
mm  downstream  of  the  injection  center.  The  images  were  taken  through  the  ceiling  of  the 
test-section,  so  a  two-dimensional  projection  of  the  curved  surface  (missile  fuselage)  was 
imaged  onto  the  CCD  camera. 

The  general  flowfield  features,  such  as  boundary  layer  separation  point, 
interaction  shock,  horseshoe  vortex,  and  the  low-pressure  wake  region,  are  labeled  in  Fig 
37.  Surface  pressure  contours  for  the  model  MIA  at  MPR  0.006  and  0.005  are  shown  in 
the  figure.  Two  high-pressure  regions  upstream  of  the  jet  were  observed,  which  is 
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consistent  with  Gruber  et  al  (1997).  The  compression  due  to  the  interaction  shock  and  the 
horseshoe  vortex  created  these  high-pressure  areas.  Upstream  of  the  interaction  shock, 
boundary  layer  separation  occurred,  and  thus  a  re-circulation  region  was  created,  which 
can  be  seen  as  the  first  thick  semicircular  high  pressure  region  upstream  of  the  jet.  For 
model  MIA,  the  surface  pressure  upstream  of  the  boundary  layer  separation  point  was 
measured  higher  (13.92  kPa)  compared  to  the  other  models,  which  was  due  to  the  higher 
tunnel  stagnation  pressure  (724  kPa)  for  this  case. 

Figures  38  -  39  show  the  surface  pressure  contours  for  models  M3A-LP  and 
M3A-HP,  respectively.  In  each  figure,  the  first  image  corresponds  to  the  low  MPR  case 
and  the  second  image  is  for  the  high  MPR  injection.  The  pressure  contours  for  all  the 
models  are  qualitatively  similar.  However,  the  pressure  in  the  separation  region  was 
observed  to  increase  with  increasing  MPR  (injection  pressure),  and  the  boundary  layer 
separation  point  moved  further  upstream  with  increased  MPR.  The  size  of  the  wake 
region  and  the  decrease  in  the  wake  pressure  both  increased  with  increasing  MPR. 

Figure  40  shows  pressure  contours  for  model  M3H  with  helium  injection. 
Because  of  the  potential  inaccuracies  due  to  the  presence  of  helium  on  the  oxygen 
sensitive  probe  molecules,  the  pressure  results  for  model  M3H  are  considered  qualitative. 
The  jet  MPR  and  jet  total  pressures  used  for  this  model  are  identical  with  model  M3A- 
LP.  As  discussed  earlier,  helium  as  a  lighter  gas  was  expected  to  simulate  higher 
temperature  air  injection.  It  can  be  gleaned  from  the  PSP  images  that,  the  boundary  layer 
separation  occurred  at  almost  the  same  location  in  absolute  distance  as  seen  with  the 
model  M3A-LP,  but  the  jet  exit  diameter  for  model  M3H  was  about  16%  lower  than 
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model  M3A-LP.  The  high-pressure  region  near  the  interaction  shock  is  much  thinner  in 
model  M3H.  A  larger  low-pressure  wake  region  is  also  observed  in  model  M3H. 

Plotted  in  Figs.  41  (a  -  d)  are  the  centerline  axial  pressure  profiles.  The  figures 
compare  the  effect  of  MPR  on  the  pressure  distribution  for  the  same  injector  geometry, 
same  injectant,  and  the  same  injection  Mach  number.  Two  high-pressure  peaks  upstream 
of  the  jet  were  observed,  as  discussed  earlier  these  peaks  resulted  from  the  compression 
due  to  the  interaction  shock  and  the  horseshoe  vortex.  Previous  studies  [Champigny  and 
Lacau  (1994),  Graham  et  al  (1999),  Hsieh  (1995),  Graham  et  al  (2000)],  both 
experimental  and  numerical,  demonstrated  similar  trends  in  the  longitudinal  surface 
pressure  distribution.  The  data  plotted  in  this  format  provides  a  clearer  documentation  of 
the  trends  listed  in  the  previous  paragraph.  Also  from  these  data,  the  recompression  shock 
appeared  to  move  downstream  for  higher  values  of  MPR. 

Figure  42  compares  the  effect  of  jet-to-freestream  pressure  ratio  on  the  surface 
pressure  for  the  models  M3A-LP  at  MPR  0.06  and  M3A-HP  at  MPR  0.05.  Due  to 
significantly  different  injector  geometry,  this  comparison  was  not  quite  justifiable. 
However,  the  shapes  of  the  curves  were  qualitatively  similar. 

The  jet  interaction  performance  is  generally  expressed  in  terms  of  the  force 
amplification  factor  [Champigny  and  Lacau  (1994)],  Kf.  The  amplification  factor  is 
defined  by  the  ratio  of  the  total  reaction  force  (F)  and  the  total  jet  delivered  force  (Fj)  at 
the  ffeestream  pressure.  For  the  present  study,  the  following  relations  were  used 


F  =  mjvJ  +(pej -pe)-Ae  +  pmg -(Alolal~Ae) 

(5.1) 

Fj=mJVj-pmAe 

(5.2) 
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(5.3) 


riijVj  is  the  jet  momentum;  pej  is  the  jet  exit  pressure;  pe  is  the  average  pressure  around 
the  circumference  of  the  exit  port;  pavg  is  the  average  pressure  on  the  entire  painted 
image;  Alolal  is  the  total  painted  area  that  was  used  to  plot  the  surface  pressure  contour, 
and  pa  is  the  freestream  pressure. 

The  jet  thrust  force  coefficient  (Cf),  defined  here  as  the  force  in  Eqn.  5.1 
normalized  by  the  freestream  dynamic  pressure  force  (pjiJ Ab  / 2  =  ypa>MlAb  / 2)  is 
plotted  in  Fig.  43  against  the  jet  MPR.  An  almost  linear  increase  in  the  jet  thrust  force 
with  increasing  of  MPR  was  evident  from  these  data,  where  CF  ~  0. 7 6 MPR  .  The  helium 

data  are  inaccurate,  and  hence  omitted;  the  presence  of  helium  results  in  a  lower 
concentration  of  oxygen  for  a  given  pressure,  which  results  in  an  erroneous  pressure 
reading.  The  amplification  factor  (KF)  was  computed  as  listed  in  Eqns  5. 1-5. 3.  The 
resulting  data  had  hyperbolic  shape.  Hence,  1/KF  is  plotted  in  Fig.  44,  and  the  trend  is 
very  nearly  linear,  where  KF  «  0.29  /  MPR . 

5.2.3  Particle  Image  Velocimetry  (PIV) 

PIV  measurements  were  performed  for  the  models  MIA  (MPR  =  0.006)  and 
M3A-LP  (MPR  =  0.06  and  0.03).  Seeder  pressure  limitations  prohibited  measurements 
for  the  underexpanded  cases.  The  PIV  measurements  were  obtained  in  x-y  planes  at  three 
spanwise  locations  (z  =  0  mm,  1.59  mm,  and  3.18  mm,  respectively).  The  mean  flow 
measurements  include  the  mean  velocity,  vorticity  and  strain  rates  (see  Chapter  4).  The 
data  were  nondimensionalized  with  the  missile  diameter  and  the  freestream  velocity. 
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The  total  number  of  images  (i.e.,  samples)  used  to  obtain  the  statistical  averages 
for  model  MIA  was  nominally  1200  (varied  from  1173-1204  across  the  three  on 
measurement  planes).  For  the  model  M3A-LP,  the  number  of  images  ranged  from 
approximately  500-800.  This  lower  number  of  samples  was  due  to  difficulties  associated 
with  pressurizing  the  seeder  for  the  Mach  3  injection.  The  contours  were  plotted  on  5120 
(80  x  64)  grid  points,  where  each  cell  was  0.85  mm  x  0.85  mm  in  size.  The  average  at 
each  point  in  the  domain  was  computed  with  at  least  70%  of  the  total  number  of  valid 
vectors. 

The  mean  velocity  contours  for  the  three  measurement  planes  are  summarized  in 
Figs.  45-50.  The  freestream  flow  is  from  left  to  right.  The  jet  flow  is  vertically  upward 
and  the  injector  center  is  located  approximately  at  x  =  0.0,  y  =  0.0,  and  z  =  0.0.  For  model 
MIA  (Figs.  45  and  46),  the  maximum  average  velocity  ratio  u  /ua  in  the  plume  was 
measured  at  approximately  0.97.  The  maximum  vertical  velocity  ratio  v  /  um  in  the  jet 
plume  was  measured  at  about  0.33  near  the  jet  exit.  Thus,  the  sonic  velocity  at  the  jet  exit 
was  confirmed.  For  model  M3A-LP  at  MPR  0.06  (Figs.  47  and  48),  the  maximum 
average  velocity  ratio  in  the  freestream  direction  was  again  about  0.97.  The  vertical 
velocity,  normalized  by  the  freestream  velocity  was  measured  as  about  0.32  after  the 
terminating  oblique  shock.  The  maximum  velocity  ratio  of  0.97  -  0.99  was  obtained  for 
the  model  M3A-LP  at  MPR  0.03  (Figs.  49  and  50)  in  the  freestream  direction. 

Focusing  on  the  mean  axial  velocity  contour  for  model  MIA  [i.e..  Fig.  45(a)],  a 
relatively  large  region  of  high  velocity  flow  was  seen  to  emanate  fromx  «  15  mm  and  y  « 
6.  The  lower  boundary  of  this  region  extends  to  the  contour  maximum  (x  «  53  mm,  y  « 
10),  and  upper  bound  reaches  the  plume  edge  at  x  «  25.  Similar  regions  with  slightly 
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different  boundaries  were  seen  for  the  two  off  centerline  planes  [Figs.  45(b)  and  (c).  The 
flow  accelerated  forward  of  this  region.  Below  this  region  (between  the  plume  core  and 
the  wake),  a  strong  velocity  gradient  was  present.  The  vertical  velocity  component  (Fig. 
46)  demonstrated  was  more  interesting.  First,  along  the  centerline  [Fig.  46(a)],  the 
vertical  velocity  was  positive  coming  out  of  the  injector  and  remained  positive  along  the 
upper  edge  until  x  «  40,  y  ~  15.  The  lower  boundary  (almost  a  line  between  extending 
from  x  «  4,  y  «  3  to  x  «  40,  y  «  1 5)  of  the  positive  vertical  velocity  is  clearly  indicated  on 
the  contour.  A  trapezoidal  region  of  negative  velocity  existed  below  the  upper  region  of 
positive  velocity;  this  region  approximately  coincides  with  the  recompression  shocks 
seen  in  the  shadowgraph.  Lastly,  below  this  negative  (downward)  velocity  region  is  a 
lower,  aft  triangular  region  of  upward  velocity.  Moving  outward  in  the  spanwise  direction 
[i.e.,  moving  from  Fig.  46(a)  to  (c)],  a  large  negative  flow  region  forms  along  the  missile 
fuselage  (y  «  0),  which  indicates  that  flow  moved  down  and  around  the  fuselage.  Models 
M3A-LP  (MPR  =  0.06  and  0.03)  displayed  similar  trends.  The  main  difference  was  the 
clearly  defined  line  of  downward  velocity  (e.g.,  the  region  centered  on  the  line  between  x 
«  25,  y  «  10  andx  «  53, y  «  18)  that  coincided  with  the  recompression  shocks. 

Mean  vorticity  contours,  normalized  by  D/uoo,  are  shown  in  Figs  51-53.  For  model 

MIA  (Fig.  51),  slightly  negative/clockwise  rotation  ( - 0.7)  about  the  z-axis  was 

observed  along  the  upper  edge  of  the  plume,  and  anti-clockwise  rotation  of 
approximately  5-10  was  seen  along  the  bottom  edge  of  the  plume.  Both  clockwise  and 
anti-clockwise  rotation  increased  at  the  off-centerline  planes.  However,  the  highest 
rotation  was  observed  near  the  trailing  edge  of  the  barrel  shock,  where  the  plume  turned 
toward  the  freestream  direction.  For  model  M3A-LP  at  MPR  0.06  (Fig.  52),  a  slight 
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rotation  (<1.0)  was  observed  along  the  plume  core.  The  highest  rotation  was  near  the 
oblique  shock  where  the  plume  suddenly  turned  toward  the  freestream  direction.  For 
model  M3A-LP  at  MPR  0.03  (Fig,  53),  the  highest  rotation  was  near  the  bottom  of  the 
plume.  Slight  (-0.2  to  -0.7)  clockwise  rotation  was  observed  near  the  upper  edge  of  the 
plume;  near  the  jet  exit  higher  clockwise  rotation  (-1.7  to  -0.5)  existed. 

The  mean  shear  strain  rate  plots,  normalized  by  D/uco,  are  shown  in  Figs  54-56. 
Common  attributes  for  all  cases  were  (1)  the  high  positive  shear  strain  rate  just  above  the 
jet  exit,  (2)  the  high  positive  shear  strain  near  the  surface  downstream  of  the  injector,  (3) 
negative  strain  along  the  upper  boundary  of  the  plume,  (4)  the  relatively  low  magnitude 
shear  in  the  upper  half  of  the  plume  core,  and  (5)  relatively  high  positive  strain 
downstream  of  the  recompression  shock.  These  strain  rate  observation  tie  closely  to  the 
turbulence  discussion  below. 

5.3  Turbulent  Flow  Measurements  (PIV) 

Turbulent  statistical  flow  properties  were  obtained  from  the  PIV  data  described  in 
Section  5.2.3.  The  statistical  properties  include  the  velocity  component  turbulence 
intensities  and  the  turbulent  shear  stresses.  The  data  were  normalized  by  the  freestream 
velocity. 

5.3.1  Turbulence  intensities 

The  turbulence  intensity  contours  for  model  MIA  (MPR  0.006)  are  given  in  Figs. 
57  and  58.  Focusing  first  on  the  axial  component  along  the  tunnel  centerline  [Fig.  57(a)], 
the  jet  flow  produced  high  levels  of  turbulence  just  above  the  jet  exit  port.  Progressing 
downstream,  the  turbulence  levels  form  three  distinct  regions.  Region  1  corresponds  to 
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the  upper  edge  of  the  plume,  were  the  turbulence  levels  were  nominally  5.0%.  Region  2 
was  defined  here  as  the  high  turbulence  region  along  the  lower  edge  of  the  plume  that 
appears  to  be  connected  to  the  high  level  turbulence  at  the  jet  exit.  In  this  region,  the 
turbulence  levels  were  approximately  15.0-20.0%.  The  third  region  was  comprised  of  the 
plume  core,  where  in  the  upper  half,  the  turbulence  levels  were  very  low  (1. 0-3.0%),  and 
moving  down  towards  the  missile  fuselage  (y  =  0),  the  level  increased  to  those  in  region 
2.  Referring  back  to  Fig.  54,  these  regions  coincided  with  regions  of  high  and  low  strain 
rate,  which  is  indicative  of  a  strain  rate  turbulence  level  causal  relationship.  Specifically, 
turbulence  was  produced  in  region  1  due  to  the  shear  layer  between  the  plume  and  the 
tunnel  crossflow.  The  low  levels  in  region  3  were  the  result  of  low  strain  (Fig.  54)  and 
potentially  turbulence  stabilization  due  to  the  rotation  associated  with  the  counter  rotating 
vortex  pair.  In  region  2,  the  strain  rates  were  the  largest,  and  thus,  the  turbulence 
production  was  expected  to  be  the  largest.  In  addition,  region  2  corresponds  to  the  low- 
pressure  wake,  which  ingested  highly  turbulent  low  momentum  boundary  layer  fluid. 

The  trends  described  for  Fig.  57(a)  prevailed  for  the  two  off-centerline  contours 
[Figs.  57(b)  and  (c)]  as  well.  However,  the  levels  in  the  three  regions  varied.  At  z—  1.59 
mm,  the  average  turbulence  was  about  5  percent  at  the  upper  edge  and  about  30.0%  along 
the  bottom  edge.  At  z  =  3.18  mm,  the  region  1  levels  were  about  10.0-15.0%,  the  region 
two  levels  were  nominally  30.0%  along  the  bottom  edge. 

The  contour  plot  of  the  vertical  component  of  the  turbulence  is  presented  in  Figs 
58  (a  -  c).  The  overall  structure  was  found  to  be  similar  to  that  in  Fig.  57.  However,  the 
transverse  levels  were  in  general  50-70%  lower  than  the  axial  levels. 
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The  turbulence  velocity  contours  for  the  Mach  3  model  M3A-LP  (MPR  0.06  and 
0.03)  is  presented  in  Figures  59  -  62.  Again,  the  overall  structure  was  similar  to  that 
described  in  regards  to  Figs.  57-58,  where  the  high  concentration  of  turbulence  at  the  jet 
exit  and  the  three  regions  were  observed. 

5.3.2  Turbulent  Stresses 

The  turbulent  stress  contours  (normalized  by  the  freestream  velocity  squared)  for 
the  sonic  injection  (model  MIA,  MPR  0.006)  are  plotted  in  Figs  63  -65.  The  normalized 

Reynolds  shear  stress  (u'v' /u2m)  contour  [Fig.  63(a)]  showed  the  same  overall  structure 

as  described  in  regards  to  Fig.  57.  However,  here  the  signs  were  significant.  Just  above 
the  injector  port  the  turbulent  shear  stresses  was  positive  at  the  peak  positive  values.  In 
region  1,  the  shear  stresses  start  out  positive,  but  progressing  downstream,  the  values 
became  negative  [this  trend  was  most  clearly  seen  in  Fig.  63(c),  however  it  was  present 
for  all  three  locations].  Region  2  contained  highly  negative  shear  stress  levels.  The 
highest  negative  values  appeared  in  downstream  region  (x  >  30  mm);  the  recompression 
shock  occurs  in  this  region.  The  second  measurement  plane  [Fig.  63  (b)]  shows  similar 
structure  with  lower  peak  intensities.  Also,  the  peak  negative  intensity  in  region  2  was 
further  upstream  than  along  the  centerline.  At  the  third  measurement  plane  [Fig.  63(c)], 
the  peak  negative  shear  stresses  in  region  2  were  further  upstream  than  for  the  second 
plane.  Hence,  the  peak  shear  levels  in  region  2  progressed  downstream  as  the  flow 
wrapped  around  the  missile.  Again,  a  strong  qualitative  correlation  between  the  shear 
strain  and  the  turbulent  shear  stress  was  observed.  The  axial  and  transverse  shear  stress 
contours  (Figs.  64  and  65)  mimic  those  for  the  turbulence  intensity  plots. 
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The  normalized  turbulent  shear  stresses  for  Mach  3  model  M3A-LP  (MPR  =  0.06 
and  0.03)  are  contoured  in  Fig.  66-71.  Overall,  the  turbulent  shear  stress  structure  for  the 
Mach  3  was  similar  to  the  sonic  model.  Aside  from  magnitude  differences,  the  only 
discrepancy  was  found  at  the  centerline  plane  for  the  model  M3A-LP  at  MPR  0.06.  The 
shear  stress  structure  shows  complicated  structure  with  both  positive  and  negative  high- 
intensities  along  the  plume  edges  near  the  surface  (region  2).  This  was  most  likely  the 
result  of  only  using  about  500  samples  to  achieve  these  second  order  statistics. 

5.3.3  Turbulent  Kinetic  Energy 

A  2-D  turbulent  kinetic  energy  was  defined  here  as  (u'2  +v'2 )/  2uJ .  The 

resultant  plots  for  Models  MIA  (MPR  =  0.006),  M3A-LP  (MPR  =  0.06)  and  M3A-LP 
(MPR  =  0.03)  are  shown  in  Figs.  72  -  74.  Because  these  TKE  was  constructed  from  the 
normal  stresses,  the  trends  in  the  TKE  echo  those  previously  discussed  in  regards  to  the 
turbulence  intensities  and  normal  stresses. 
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6.  CONCLUSIONS  AND  PUBLICATIONS 


The  development  flight  control  methods  for  Army  supersonic  missiles  at  low 
altitudes  using  low  blowing  ratio  reaction  jets  is  complicated  by  (1)  nonlinear, 
counterintuitive  surface  interactions,  (2)  hot  gas  effects  and  (3)  poor  predictability  using 
modem  computation  fluid  dynamic  codes,  where  turbulence  modeling  is  a  key  factor 
limiting  the  accuracy.  The  objectives  of  this  research,  which  were  directly  tied  to  the 
Army  low  blowing  ratio  reaction  control  jet  requirements,  were  documentation  and 
improved  understanding  of  (1)  the  surface  flow  interactions  with  hot  gas  effects  and  (2) 
the  turbulent  flow  structure  in  the  plume.  To  meet  these  objectives,  an  experimental 
analysis  of  the  mean  and  turbulent  flow  properties  for  sonic  and  Mach  3  injection  into  a 
Mach  3  freestream  was  accomplished.  The  range  of  conditions  included:  MPR  =  0.005- 
0.09,  Pj/pK=  0.6-19.0,  Tj/Tk  =  1.0-13.6  (helium  was  used  to  simulate  the  highest 

temperature),  y  =  1.4  and  1.67,  and  peb/pm  =  0.19-6.15  (overexpanded,  perfectly 
expanded  and  underexpanded  jets  were  tested).  Numerical  simulations  were  not 
accomplished  in  the  present  effort.  However,  a  full  characterization  of  the  high-speed  jet 
interaction  flowfield  over  a  blended  body  missile  fuselage  (sharp-coned  cylinder)  was 
experimentally  achieved,  where  specific  pragmatic  and  scientific  flowfield  conclusions 
were  drawn.  The  conclusions,  which  are  separated  into  three  categories  based  on 
measurement  technique  (flow  visualization,  mean  flow,  and  turbulent  flow),  are 


summarized  below. 


6.1  Flow  Visualization 


Quantitative  and  qualitative  information  concerning  the  flow  topology,  e.g.,  the 
locations  of  the  interaction  shock,  horseshoe  vortex,  recompression  shock,  and  wake 
region,  were  discerned  from  the  flowfield  visualizations.  From  these  data,  the  following 
observations  were  deemed  important.  First,  the  surface  oil  flow  visualization  provided 
qualitative  insight  into  the  surface  topology  and  quantitative  measurement  of  the 
upstream  separation  distance,  which  have  important  implications  relative  to  the 
amplification  factor  for  flight  missiles.  From  the  surface  oil  flow  visualizations,  the 
following  conclusions  were  drawn: 

•  Injection  temperature,  for  a  given  MPR,  had  a  small  but  discernible  effect  on 
the  boundary  layer  separation  distance,  where  for  a  given  MPR,  the  separation 
distance  increased  with  temperature  ratio  at  a  rate  of  approximately  0.013. 
This  increase  was  smaller  than  expected  from  previous  overexpanded  studies, 
however  the  present  high  temperature  cases  were  all  overexpanded,  and  hence 
the  upstream  influence  was  smaller. 

•  The  boundary  layer  separation  increased  significantly  with  MPR,  however 
MPR  was  not  a  sufficient  correlating  parameter.  For  example,  for  a  given 
MPR,  the  separation  distance  increased  when  the  jet  condition  transitioned 
from  overexpanded  to  underexpanded. 

•  The  influence  of  the  jet  plume  on  the  surface  pressure  distribution  in  the  wake 
region  of  the  flow  was  found  to  depend  on  MPR,  where  the  wake  region 
developed  further  downstream  in  terms  of  distance  nondimensionalized  by  the 
jet  exit  diameter. 
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Second,  the  shadowgraph  and  schleiren  photographs  provided  documentation  of  the 
shock  structure  along  the  flow  centerline.  The  following  observation  were  drawn  from 
the  shadowgraph  and  schlieren  images: 

•  The  effective  backpressure  concept  was  an  accurate  indicator  of  the  jet  flow 
condition  (overexpanded  versus  underexpanded),  where,  for  example,  a 
terminating  Mach  cone  was  observed  for  all  of  the  cases  that  were  deemed 
overexpanded  using  this  criterion. 

•  The  jet  interaction  bow-shock  strength  increased  with  increased  MPR. 

•  For  a  given  MPR,  overexpanded  jets  produced  weaker  shocks  than 
underexpanded  jets,  where  the  overexpanded  shock  structure  turned  the  flow 
downstream,  thus  lessoning  the  flow  disturbance. 

Third,  the  Mie-scattering  flow  visualizations  provided  qualitative  insight  into  the 
turbulent  structure  of  the  plume,  the  location  of  the  plume  surface  interaction  and  the 
penetration  of  the  plume  into  the  crossflow.  The  following  observations  were  deemed 
important: 

•  Large-scale  turbulent  structures  were  present  in  the  plumes. 

•  The  surface  interaction  location  moved  closer  to  the  jet  exit  with  decreasing 
MPR. 

•  The  penetration  height  increased  with  increasing  MPR. 
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6.2  Mean  Flow 


Detailed  mean  flow  documentation  was  achieved  with  pressure  sensitive  paint, 
the  five-hole  pressure  probe  and  particle  image  velocimetry.  First,  The  PSP  images 
provided  both  qualitative  flow  visualization  and  quantitative  mean  surface  static  pressure 
measurements.  From  the  PSP  data  the  following  observations  were  drawn: 

•  The  upstream  peak  pressures  associated  with  lambda  shock  and  the  horseshoe 
vortex  increased  with  increasing  MPR. 

•  The  wake  pressures  decreased  with  increasing  MPR  (or  injection  pressure). 

•  The  interaction  force  increased  with  almost  linearly  increasing  MPR  (i.e., 
CF  &0.76- MPR ). 

•  The  amplification  factor  was  inversely  proportional  to  MPR  (i.e., 
Kf»  0.29 /MPR). 

Second,  the  five-hole  probe  provided  the  detailed  contours  of  Pitot  pressure,  Mach 
number,  total  pressure,  and  static  pressure.  From  the  five-hole  probe,  data  the  following 
additional  conclusions  were  drawn 

•  The  mean  flow  features  (plume,  wake  region,  interaction  shock,  counter 
rotating  vortex  pair  and  the  horseshoe  vortex)  were  documented. 

•  The  counter-rotating-vortex  pair  became  more  prominent  with  increasing 
MPR. 
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Third,  the  PIV  provided  the  contours  of  the  mean  velocity,  z-component  of  vorticity,  and 
mean  x-y  shear  strain  rate  at  three  z-locations  for  models  MIA  (MPR  =  0.006),  M3A-LP 
(MPR  =  0.03  and  0.06).  From  the  PIV  data,  the  following  observations  were  drawn: 

•  The  mean  axial  and  transverse  velocity  components  across  the  plume  were 
influenced  by  both  the  plume  vorticity  and  surface  curvature,  where  in  the 
upper  region  of  the  plume  where  the  counter  rotating  vortex  pair  resided,  the 
axial  velocity  peaked  a  value  nearly  that  of  the  freestream,  and  the  data 
showed  that  jet  fluid  moved  down  and  around  the  missile  fuselage. 

•  Clockwise  vorticity  (z-component)  was  observed  along  the  upper  edge  of  the 
plume,  while  anti-clockwise  vorticity  was  observed  near  the  bottom  edge  of 
the  plume. 

•  The  shear  strain  rate  fields  were  complicated,  where  a  region  of  high  positive 
shear  strain  rate  was  seen  in  the  jet  plume  just  above  the  injector  port,  where 
the  plume  first  interacted  with  the  high  momentum  crossflow.  A  layer  of  high 
shear  strain  rate  existed  along  the  missile  fuselage.  The  shear  strain  had 
relatively  low  magnitudes  in  the  upper  half  of  the  plume  core,  and  it  increased 
across  the  recompression  shock. 

6.3  Turbulent  Flow 

Turbulent  flowfield  was  documented  with  PIV  for  models  MIA  (MPR  =  0.006), 
and  model  M3A-LP  (MPR  =  0.03  and  0.06).  The  measurements  included  velocity 
component  turbulent  intensities,  turbulent  stresses,  and  turbulent  kinetic  energy.  From 
these  data,  the  following  conclusions  were  made: 
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•  In  general,  progressing  downstream,  the  turbulence  levels  formed  three 
distinct  regions  across  the  plume. 

o  Region  1  was  defined  as  the  upper  region  of  the  plume  where  the 
turbulence  levels  were  moderate  (turbulence  intensities  ~  5.0-10.0%). 
In  this  region,  the  Reynolds  shear  stress  started  out  positive  and 
transitioned  to  negative  progressing  downstream, 
o  Region  2  was  defined  as  the  high  turbulent  level  (turbulence  intensities 
~  15.0-30.0%)  region  located  along  the  lower  portion  of  the  plume 
near  the  missile  fuselage.  The  Reynolds  shear  stresses  had  large 
magnitude  negative  values  in  this  region, 
o  Region  3  was  defined  as  the  relatively  low  turbulence  region  along  the 
plume  core,  where  in  the  upper  half,  the  turbulence  levels  were  very 
low  (turbulence  intensities  ~  1. 0-3.0%),  and  moving  down  towards  the 
missile  fuselage,  the  levels  increased  to  those  in  region  2. 

•  The  moderately  high  turbulence  levels  in  region  1  were  due  production  that 
resulted  from  the  strain  rates  across  the  shear  layer  between  the  plume  and  the 
tunnel  crossflow. 

•  In  region  2,  the  strain  rates  were  the  largest,  and  thus,  the  turbulence 
production  was  expected  to  be  the  largest.  In  addition,  region  2  corresponded 
to  the  low-pressure  wake,  which  ingested  highly  turbulent  low  momentum 
boundary  layer  fluid. 
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•  The  low  levels  in  region  3  were  the  result  of  low  strain  rates  and  turbulence 
stabilization  due  to  the  rotation  associated  with  the  counter  rotating  vortex 
pair. 

•  The  present  data  suggest  a  strong  qualitative  correlation  between  the  shear 
strain  rate  field  and  the  turbulence  across  the  plume.  This  finding  supports  the 
usage  of  the  Boussinesq  turbulence  model  approximation. 

6.4  Publications 

The  following  two  conference  papers  were  generated  from  a  subsection  of  the 
results  described  above: 

•  Mahmud,  Z.  and  Bowersox,  R.,  “Supersonic  Missile  Body  Jet  Interaction 
Flowfields  at  Low  Momentum-Ratio-Parameter,”  to  be  presented  at  the  41st 
AIAA  Aerospace  Sciences  Meeting,  Reno  NV,  Jan  2003. 

•  Mahmud,  Z.,  and  Bowersox,  R.,  "Experimental  Investigation  of  Low  Blowing 
Ratio  Blended  Body  Supersonic  Injection  Flowfields,"  AIAA  2001-0886,  39th 
AIAA  Aerospace  Sciences  Meeting,  Reno  NV,  Jan  2001. 
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APPENDIX  A  FREESTREAM  UNIFORMITY 


The  five-hole  probe  was  used  to  measure  the  freestream  Mach  number  contour  (Fig.  75). 
The  average  Mach  number  was  2.97±2.0%.  The  freestream  turbulence  level  was 
measured  as  3  -  4%  using  the  PIV  technique.  The  freestream  was  seeded  through  a  hole 
in  the  settling  chamber.  The  seeder  pressure  was  about  20  -  30%  higher  than  the  settling 
chamber  pressure.  Hence,  the  accuracy  of  this  measurement  is  suspect. 
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APPENDIX  B  UNCERTAINTY  ANALYSIS 


B.l  Five-Hole  Probe 

Pressure  transducer  uncertainties  include  manufacturer  uncertainty  (linearity  & 
hysteresis,  and  repeatability,  calibration  uncertainty,  and  resolution  uncertainty  of  the 
data  acquisition  card.  The  manufacturer  uncertainties,  as  tabulated  in  Table  15,  are 
converted  into  pressure  values.  Calibration  uncertainties  were  computed  from  the 
differences  between  the  actual  pressure  reading  and  the  estimated  pressure  values  using 
the  calibration  equations.  The  resolution  uncertainty  of  the  acquisition  card  of  0.0024 
volts  was  converted  into  the  corresponding  pressure  value  using  the  transducer 
calibration  equations.  The  transducer  calibration  equations  are  described  below. 

pt2  (psiA)  =  25.029*  Volts  -  6.2268 
PcI  (psiA)  =  3. 7421  *  Volts  -  0.5223 
pc2  (psiA)  =  3. 7724  *  Volts  -  0.9121 
pC3  (psiA)  =  3.7663  *  Volts  -  0.6494 
pc4  (psiA)  =  3.7681  *  Volts  -  0.4492 
Ptoo  (psiA)  =  19.19  *  Volts  -  0.0866  +  patm 
Pec  (psiA)  =  2.9551  *  Volts  +  0.3303 

Detailed  transducer  uncertainties  are  listed  Table  15.  The  Thermocouple 
uncertainties  are  listed  in  Table  16. 
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The  total  uncertainty  in  Mach  number  measurement  includes  the  uncertainties  due 
to  the  pressure  transducer,  5-hole  probe  Mach  number  calinration,  and  probe  inclination. 
The  pressure  transducer  uncertainties  are  described  above.  The  Mach  number  uncertainty 
of  ±  0.03  is  documented  for  the  probe  misalignment  of  ±  6  degree  [Volluz  (1961)].  The 
overall  L2-norm  of  the  uncertainty  was  computed  following  Holman  (1978);  that  is 
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where,  M,  =  Mach  number  from  the  Taylor-Mccoll  conical  from  theory  approximation; 
coM.  =  Uncertainty  for  M' (=  1/M, );  o/M  =  Uncertainty  for  M, ;  Mrcal  =  Mach  number 
based  on  experimental  re-calibration;  corcal  =  Uncertainty  for  Mrcal;  coangle  =  Uncertainty 
of  Mach  number  for  probe  inclination  angle;  cocaI  =  Uncertainty  for  the  curve-fit  solution 
of  the  conical  flow  theory,  and  coM  (  =  Total  uncertainty  of  Mach  number  from  the  5. 
The  inviscid  solution  was  curve  fit  by  Eqn  (4.1),  the  hence  the  derivatives  were  given  by 
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The  Pitot  pressure  uncertainty  includes  both  probe  angle  effects  and  transducer 
calibration.  The  static  and  total  pressure  uncertainties  were  computed  following  the  above 
analysis  using  the  compressible  flow  relation.  The  resulting  uncertainties  are  listed  in 
Table  17. 

B.2  LaVision  PIV 

The  PIV  uncertainty  stemmed  from  the  assumption  that  the  particle  position 
uncertainty  was  0.3  pixels  [Gogineni  et  al  (1998)].  With  this,  the  uncertainties  for  the 
measured  quantities  were  estimated  as  listed  in  Table  17. 

B.3  Pressure  Sensitive  Paint 

The  pressure  sensitive  paint  uncertainty  estimate  included  system  (calibration  and 
camera)  as  provided  by  the  manufacturer  (2.0%  of  the  measured  value),  and  surface 
temperature.  The  surface  temperature  was  known  to  approximately  ±5K,  and  the 
temperature  sensitivity  was  0.6%/K.  Hence,  the  overall  L2  norm  was  computed  as  3.6% 
of  the  measured  value. 

B.4  Separation  Distance 

The  separation  distance  was  measured  with  a  dial  caliper  that  has  a  resolution  of 
25.4  pm.  This  resolution  was  taken  as  the  measurement  uncertainty. 
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TABLES 


_ Table  1  Range  of  Jet  Parameters 

Jet  Mach,  M,  MPR _ Pj/P^  T/  To, 

1.0  -3.0  0.005-0.09  0.6-  18.8  1.0-  13.6  1.4,1.67 
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Table  2  Freestream  Conditions 


Mach 

Velocity,  m/s 

Temperature,  K 

Pressure,  KPa 

Re/m 

3.0 

616 

105 

17.3 

48.5 

3.0% 
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Table  3  Tunnel  Operating  Conditions 


Ptoo,  Mpa 

Ttoo,  K 

0.6067  ±  0.0207 

285  ±  2.2 
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Table  4  Diffuser  Dimensions 


HI  Mach 

A 

Hl-A 

<t> 

LI 

Ll/W 

n 

L2  H2 

L2/W 

a 

H 

3  4.0 

2.00 

1 

9.50 

6 

3.006 

8.98 

19 

5 

9.52 

14.07 

4.41 

1.48 

3.0 

2.14 

0.86 

8.15 

6 

2.802 

8.56 

19 

5 

8.87 

12.72 

3.98 

1.34 

2.0 

2.45 

0.55 

5.26 

6 

2.451 

7.66 

19 

5 

7.76 

9.83 

3.08 

1.03 

1.5 

2.73 

0.27 

2.61 

6 

2.2 

6.82 

19 

5 

6.97 

7.18 

2.25 

0.75 

(see  figure  3.5  for  clarity) 


I 
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Model 

Exit 

Mach 

Throat  diameter,  Dj, 
mm  (in) 

Exit  diameter,  De, 

mm  (in) 

De/Dj 

M3A-HP 

3.0 

1.191  (3/64) 

2.453 

2.06 

M3A-LP 

3.0 

3.969  (5/32) 

8.176 

2.06 

M3H 

3.0 

3.969  (5/32) 

6.863 

1.73 

MIA 

1.0 

4.763  (3/16) 

4.763 

1 

M1A-L 

1.0 

9.525  (3/8) 

9.525 

1 
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Table  6  Injector  Operating  Conditions 


Model 

M} 

Injectant 

Ptj 

MPa 

Pj/Poo 

Pj/Peb 

Yj 

YjM// 

Yoo  Moo 

At/Aj 

MPR 

M3A- 

HP 

3.0 

Air 

5.96 

9.38 

3.07 

1.4 

1.00 

210.2 

0.05 

12.9 

18.82 

6.15 

0.09 

M3A- 

LP 

3.0 

Air 

0.377 

0.60 

0.19 

1.4 

1.00 

18.9 

0.03 

0.694 

1.09 

0.36 

0.06 

M3H 

3.0 

He 

0.365 

0.66 

0.22 

1.67 

1.19 

26.9 

0.03 

0.670 

1.21 

0.40 

0.05 

MIA 

1.0 

Air 

0.088 

2.68 

0.88 

1.4 

0.11 

55.8 

0.005 

0.101 

3.10 

1.01 

0.006 

0.263 

8.04 

2.06 

0.015 

M1H 

1.0 

He 

0.224 

6.31 

2.06 

1.67 

0.13 

55.8 

0.015 

M1A-L 

1.0 

Air 

0.088 

2.68 

0.88 

1.4 

0.11 

13.9 

0.021 

0.101 

3.10 

1.01 

0.024 

MPR=  YjMj2PjAj  /YooMoo2PooAb;  Peb  =  0.8P2;  A  =  Air;  H  =  Helium 


0.024 


Model 

MPR 

Tu,K 

M3A-HP/LP 

0.03-0.06 

294.0 

1.0 

MIA 

0.005-0.015 

294.0 

2.3 

M1A-L 

0.021,0.024 

294.0 

2.3 

MIA 

0.006 

920.0 

7.3 

M3H 

0.006 

1250.0b 

11.9 

M1H 

0.015 

1 900.0b 

13.6 

“Pebble  bed  heated  air  injection,  bSimulated  temperature  based  on  light  gas  (helium)  injection 


1 
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_ Table  8  Flowmaster3  Camera  Specs 

Model  Flowmaster3S 

Sensor  SONY  ICX  085 


Pixel  (h  x  v) 

Pixel  Size  (h  x  v) 
Active  Area  (h  x  v) 
Fill  Factor 
Dynamic  Range 
Single  Frames  per 
Second 

Min  Interframing 
Time 


1280  x  1024  (square 
format) 

6.7  pm  x  6.7  pm 
8.6  mm  x  6.9  mm 
60%  (with  microlens) 

12  bit 

8 

300  ns 
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Wavelength  532  nm 

Gemini  PIV1 5  120mJ 

Energy  Stability  ±3 . 5  % 

Pulse  width  3-5ns 

Beam  divergence  <1  mrad 

Beam  pointing  <200  prad 

Jitter  ±0.5  ns 


_ Table  10  Transducer  Specifications 

PX139  series  PX139  series 


Excitation  Voltage 
Output: 

Linearity  and  Hysteresis: 
Repeatability: 

Proof  pressure: 

Burst  Pressure: 


5  Vdc 

0.25  to  4.25  volts 
±0.1%  FS 
±0.3%  FS 

>  3  x  FS  Pressure 

>  5  x  FS  Pressure 


PX202-100GV,  and  PX302-015AV  PX202-100GV,  and  PX302-015AV 


Excitation  Voltage 
Output: 

Linearity,  Hysteresis,  and  Repeatability 
Response  Tine 

Vibration: 

Proof  pressure: 

Burst  Pressure: 


10  Vdc 

lOOmV  ±  lmV 
±  0.25%  FS 
lmsec 

PX202  -  35g  peak  @  5  -  2000Hz 
PX302  -  15g  peak  @  10  -  2000Hz 
PX202  -  150% 

PX302  -  200% 

400% 
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Table  11  Five-Hole  Probe  Calibration 


Mach 

%error  2 

1.20 

1.30 

0.73 

1.40 

0.498 

0.37 

1.50 

0.4568 

0.47 

1.75 

0.3794 

0.62 

2.00 

0.3256 

0.20 

2.25 

0.2868 

0.07 

2.50 

0.258 

0.02 

2.75 

0.2359 

0.25 

3.00 

0.2187 

0.50 

3.25 

0.2051 

0.68 

3.50 

0.1941 

0.71 

4.00 

0.1777 

0.42 

4.50 

0.1661 

0.36 

5.00 

0.1580 

0.79 

6.00 

0.1470 

1.13 

7.00 

0.1402 

1.21 

8.00 

0.1357 

0.60 

9.00 

0.1325 

0.74 

10.00 

0.1300 

0.91 

’Theory:  Exact  for  Mach<5  and  10.0,  Hypersonic  Similarity  for  Mach=6-9 
(Difference  between  exact  and  Hypersonic  Similarity<0.5%  at  Mach=5.0  and  10.0). 

2  based  on  the  curve-fit  solution. 
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Table  12  Five-Hole  Probe  Calibration  Data 


Pc/Pt2 

Mach1 

Mach2 

Mach 

0.142 

7.04 

5.52 

5.55 

0.155 

6.45 

4.83 

4.79 

0.229 

4.37 

3.00 

3.01 

0.452 

2.21 

1.72 

1.72 

conical  flow  theory,  normal  shock  relationship, Recalibrated  curve  fit 


I 
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Table  13  Parameters  used  to  process  PIV  images  and  vectors 


Image  processing: 

_ Mode:  Cross-correlation  (1  camera) _ 

Int.  window  size  and  shift: _ 

Adaptive  multi-pass  (4  passes)  with  decreasingly  smaller  sizes  were  used 

Initial  int.  window  size:  64  x  64 

Final  int.  window  size:  32  x  32  [iteration:  3] 

Overlap:  50% 

_ Initial  shift:  0 _ 

Correlation  Function  [Initial  and  Final  passes): _ 

Standard  11*12  (via  FFT,  no  zero  padding) 

_ Deformed  int.  window:  off _ 

Multi-pass: _ 

Relative  vector  range  restriction:  reference  ±  cell-size  /  8 
AND  absolute  vector  range  restriction:  reference  ±  4  pixels 

_ Overlap  for  initial  passes:  50% _ 

VECTOR  POST  PROCESSING _ 

Median  Filter: 

_ Strongly  remove  and  iteratively  replace _ 

Remove  if  >  1.5  *  RMS  of  neighbors 
Remove  if  <  3  neighbors 

_ Insert  if  <  2  *  RMS  of  neighbors _ 

_ Remove  groups  with  less  than:  10  vectors _ 

_ Smooth:  3x3 _ 

_ Allowed  vector  range: _ 

Vx:  -400  to  700  m/s 

Vy:  -400  to  700  m/s  for  Mach  3  injection  case 
Vy:  -400  to  400  m/s  for  sonic  injection  case 


Table  14  Boundary  Layer  Separation  Distance 


Separation  Distance 

Model 

Tj/Too 

De,  mm 

MPR 

Pj/Pinf 

(A/De) 

M3A-HP 

1.0 

2.45 

0.09 

0.05 

18.82 

9.38 

2.81 

2.08 

M3A-LP 

1.0 

8.18 

0.06 

0.03 

1.09 

0.59 

1.88 

1.61 

M3H 

11.9 

6.86 

0.05 

0.03 

1.21 

0.66 

2.22 

1.75 

MIA 

2.3 

0.02 

3.10 

1.76 

4.76 

0.02 

2.68 

1.75 

MIA 

7.3 

0.02 

3.10 

1.83 

MIA 

2.3 

4.76 

0.015 

8.13 

2.03 

M1H 

13.6 

4.76 

0.015 

6.41 

2.16 
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Table  15  Pressure  Transducers  and  Uncertainties 


Transducer 

Model 

Measured 

pressure 

Manufacturer 
uncertainty 
+  Pa  (psia) 

Calibration 
uncertainty 
+  Pa  (psia) 

Resolution 
uncertainty  of 
acquisition  card 
±  Pa  (psia) 

Total 

Uncertainty 
±  Pa  (psia) 

Omega 

PX139- 

100A4V 

Pa 

2180  (0.316) 

446.8 

(0.0648) 

413.7  (0.06) 

2260  (0.328) 

Omega 
PX139- 
01  SAW 

Pci 

326.8 

(0.0474) 

92.4 

(0.0134) 

62.1  (0.009) 

344.7  (0.05) 

Omega 
PX139- 
01  SAW 

Pc2 

326.8 

(0.0474) 

139.3 

(0.0202) 

62.1  (0.009) 

358.5  (0.052) 

Omega 
PX139- 
01 5 AW 

Pci 

326.8 

(0.0474) 

122.7 

(0.0178) 

62.1  (0.009) 

351.6  (0.051) 

Omega 
PX139- 
01 5 AW 

Pc4 

326.8 

(0.0474) 

147.6 

(0.0214) 

62.1  (0.009) 

365.4  (0.053) 

Omega 

PX202- 

100GV 

Ptcc 

1720  (0.25) 

176.3 

(0.02557) 

317.16 

(0.046) 

1758  (0.255) 

Omega 

PX302- 

015AV 

Poo 

258.6 

(0.0375) 

93.08 

(0.0135) 

46.88 

(0.0068) 

275.8  (0.04) 

Table  16  High  Temperature  Thermocouple  and  Uncertainty 


Thermocouple  Model 

Measurement  Uncertainty  ±  °C  (°F) 

Omega  XTA-W5R26-  U-125-30-H-HX-6 

Ty  0.4  °C  (0.7  °F) 

78 


Table  17  Measurement  Uncertainty  Summary 

Measurement 

Uncertainty 

M  (Five-Hole  Probe) 

0.13 

Ptf/PtoO 

0.02 

Pl/Ptoo 

0.15 

Ptl/Ptoo 

0.10 

ui  /  uK  (LaVision  PIV) 

0.02 

0.02 

u'u'/u 2 

1  J  ao 

0.001 

&uj  /  dXj 

1.3 

p(PSP) 

3.6% 

A  (pm) 

25.4 
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FIGURES 


Figure  1  Low-Speed  Jet  Interaction  Flowfield 
[Margason  (1993)] 
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20 


14 


1.  Compressors 

2.  Supply  tank 

3.  Hand  operated  valve 

4.  Electric  valve 

5.  Hydraulic  Pressure  regulator 

6.  Flow  straightener 

7.  Mach  3  nozzle 

8.  Test  section 

9.  Missile  model 

10.  inserted  probe 

11.  diffuser 

12.  probe  traversing  motor 

13.  air  bottle 

14, 15, 16.  tables 

17.  on-off  switch  (electric  valve) 

18.  Data  acquisition  and  computer 
assembly 

19.  Jet 

20.  exhaust  air 

21.  Horizontal  traverse  of  the  probe 
assembly 

22.  vertical  traverse  of  the  probe 


(a)  Schematic 


(b)  Photograph 

Figure  3  Mach  3.0  Wind  Tunnel  (schematic  and  photograph) 
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Figure  4  Flow  Straightener  (schematic) 
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Figure  5  Sketch  of  the  Mach  3  Nozzle 
(Not  to  Scale) 


(a)  Schematic 


(b)  Photograph 

Figure  6  Test  section  with  a  model  (schematic  and  photograph) 
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(a)  Electric  Furnace 


(b)  Hot  Ceramic  Pebbles 


heater 


(d)  Pebble  Bed  and  Tunnel 
Figure  11  Pebble  Bed  Heater 
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(a)  PIV  Seeder  Tank 


(b)  LaVision  PIV  system 
Figure  13  LaVision  PIV  Apparatus 
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Figure  14  Five-Hole  Probe  Calibration 
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Interaction 
shock  region 


Jet 


Horseshoe 

vortex 


Figure  15  General  Flow  Features 


(a)  Model  MIA,  MPR  0.006 


(c)  Model  MIA,  MPR  0.006,  at  922  K 
Figure  16  Surface  Oil  Flow  Visualizations  for  Model  MIA 
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(a)  Model  M3A-LP,  MPR  0.03 


j  . 

am  1 

'  Y£i 

• HP  •  •-  ■  SL. 

r'  f  I 

(b)  Model  M3A-LP,  MPR  0.06 

Figure  17  Surface  Oil  Flow  Visualizations  for  Model  M3A-LP 
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(a)  Model  M3A-HP,  MPR  0.05 


(a)  Model  M3A-HP,  MPR  0.09 

Figure  18  Surface  Oil  Flow  Visualizations  for  Model  M3A-HP 
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(a)  Model  M3H,  MPR  0.03 


(b)  Model  M3H,  MPR  0.06 


(c)  MPR  0.015  (air  and  helium  injection) 

Figure  19  Surface  Oil  Flow  Visualizations  for  Model  M3H 


(i)  MIA,  air  injection 


(ii)  M1H,  Helium  injection 
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Figure  22  Change  in  Boundary  layer  separation  distance  vs.  Change  in 

Temperature 


Figure  23  Mie-scattering  image,  Model  MIA,  MPR  0.006 
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(a)  Model  MIA  at  z  =  0 


(b)  Model  MIA  at  z  =  1/16  inch 


(c)  MIA  at  z  =  2/16  inch 

Figure  24  Instantaneous  Mie-scattering  images  (Model  MIA,  MPR  0.006) 
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(a)  Average  PIV  image  at  station  z  = 


I 


(b)  Average  PIV  image  at  station  z  =  1.59  mm 


I 

i 


(c)  Average  PIV  image  at  station  z  =  3.18  mm 
Figure  25  Average  Mie-scattering  images  (Model  MIA,  MPR  0.006) 


(a)  Instantaneous  particle  image  at  station  z  = 


(b)  Instantaneous  particle  image  at  station  z  =  1.59  mm 


(c)  Instantaneous  particle  image  at  station  z  =  3.18  mm 
Figure  26  Instantaneous  Mie-scattering  images  (Model  M3A-LP,  MPR  0.06) 
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ft 


(a)  Average  PIV  image  at  station  z  =  0 


(b)  Average  PIV  image  at  station  z  =  1.59  mm 


(c)  Average  PIV  image  at  station  z  =  3.18  mm 
Figure  27  Average  Mie-scattering  images  (Model  M3A-LP,  MPR  0.06) 
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(a)  Instantaneous  Mie-scattering  at  station  z  =  0 


(b)  Instantaneous  Mie-scattering  at  station  z  =  1.59  mm 


(c)  Instantaneous  Mie-scattering  at  station  z  =  3.18  mm 
Figure  28  Instantaneous  Mie-scattering  images  (Model  M3A-LP,  MPR  0.03) 
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(a)  Average  PIV  image  at  station  z  = 


0 


i 


(a)  Model  MIA,  MPR  0.005 


(b)  Model  MIA,  MPR  0.006 


(c)  Model  M1A-L,  MPR  0.003 


(d)  Model  M1A-L,  MPR  0.006 


(e)  Model  M3A-LP,  MPR  0.03 


(f)  Model  M3A-LP,  MPR  0.06 


(a)  Model  MIA,  MPR  0.005  (b)  Model  MIA,  MPR  0.006 


(c)  Model  M3A-LP,  MPR  0.03  (d)  Model  M3A-LP,  MPR  0.06 


(e)  Model  M3A-HP,  MPR  0.05  (f)  Model  M3A-HP,  MPR  0.09 

Figure  31  Color  Schlieren  Images 


Interaction! 

Shock 


3 


MPR  0.005 


7 


(b)  Model  MIA,  MPR  0.005 
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-LP,  MPR  0.06 


22 


Freestream 


Freestream 


Separation 


Horseshoe 

Vortices 


Low  pressure 


(a)  MPR  0.006 


(b)  MPR  0.005 


Figure  37  Surface  pressure  contour  (Model  MIA) 


2  3  4  5  6 


0  12  3  4  5  6 


psiA 


psiA 


(a)  MPR  0.03  (b)  MPR  0.06 

Figure  38  Surface  pressure  contour  (Model  M3A-LP) 
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0123456  012345 


psiA 


psiA 


(a)  MPR  0.05  (b)  MPR  0.09 

Figure  39  Surface  pressure  contour  (Model  M3A-HP) 
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3  4  5  6 


0  1  2  3  4  5 


psiA 


psiA 


(a)  MPR  0.03  (b)  MPR  0.05 

Figure  40  Surface  pressure  contour  (Model  M3H) 
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Model  M3A-HP 


Figure  41  Effect  of  jet  MPR  on  the  centerline  surface  pressure 


x,  mm 


(c)  z  =  3.18  mm 

Figure  45  Mean  axial  velocity  contours,  u  /«„  (Model  MIA,  MPR  =  0.006) 
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(a)  z  =  0 


-0.08-0.07-0.05-0.04-0.02-0.01  0.01  0.02  0.04  0.05  0.07  0.08  0.10  0.11  0.13  0.14  0.16  0.17  0.19  0.20 


20  30 

x,  mm 

(b)  z  =  1.59  mm 


-0.07-0.06-0.05-0.04-0.02-0.01-0.000.01  0.02  0.03  0.05  0.06  0.07  0.08  0.09  0.10  0.12  0.13  0.14  0.15 


0  10  20  30  40  50 

x,  mm 

(c)  z  =  3.18  mm 

Figure  46  Mean  transverse  velocity  contours,  v/ux  (Model  MIA,  MPR  =  0.006) 
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x,  mm 

(b)  z  =  1.59  mm 
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0.00  0.05  0.11  0.16  0.21  0.26  0.32  0.37  0.42  0.47  0.53  0.58  0.63  0.68  0.74  0.79  0.84  0.89  0.95  1.00 


x,  mm 

(c)  z  =  3.18  mm 

Figure  47  Mean  axial  contours,  « /Ua>  (Model  M3A-LP,  MPR  =  0.06) 
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x,  mm 

(c)  z  =  3.18  mm 

Figure  48  Mean  transverse  velocity  contours,  v  /  ux  (Model  M3A-LP,  MPR  =  0.06) 
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Figure  49  Axial  Velocity  Contours,  w/wM  (Model  M3A-LP,  MPR  =  0.03) 
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Figure  51  Vorticity,  (du  / dy-dv /dx)  D/ux  (Model  MIA,  MPR  =  0.006) 
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Figure  52  Vorticity,  (du  / dy-dv  / dx)- D /u„  (Model  M3 A-LP,  MPR  =  0.06) 
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Figure  54  Shear  Strain,  (du  /dy  +  dv  / dx)-D / ux  (Model  MIA,  MPR 
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Figure  55  Shear  Strain,  (du/dy  +  dv/dx)-D/ux  (Model  M3A-LP, MPR  =  0.06) 
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Figure  56  Shear  Strain,  (dU /dy+dt /dx)  D/u„  (Model  M3A-LP,  MPR  =  0.03) 
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Figure  58  Transverse  Turbulence  Intensity, 


um  (Model  MIA,  MPR  =  0.006) 
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Figure  60  Transverse  Turbulence  Intensity,  ^/v'Vwo^  (Model  M3A-LP,  MPR  0.06) 
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Figure  61  Axial  Turbulence  Intensity,  /um  (Model  M3A-LP,  MPR  0.03) 
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Figure  62  Transverse  Turbulence  Intensity,  ^/v'^/wo:  (Model  M3A-LP,  MPR  0.03) 
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Figure  63  Reynolds  Stress,  u'v'/uj  (Model  MIA,  MPR  0.006) 
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Figure  64  Axial  Stress,  u'u'/uj  (Model  MIA,  MPR  0.006) 


154 


(a)  z  =  0 


O.OOE+OO  9.47E-04  1.89E-03  2.84E-03  3.79E-03  4.74E-03  S.68E-03 


20  30 

x,  mm 

(b)  z  =  1.59  mm 


10 
5 
0 

0  10  20  30  40  50 

x,  mm 

(c)  z  =  3.18  mm 

Figure  65  Transverse  Stress,  v'v'/uj  (Model  MIA,  MPR  0.006) 
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Figure  66  Reynolds  Stress,  u'v'/uj  (Model  M3A-LP,  MPR  0.06) 
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Figure  67  Axial  Stress,  u'u' / uj  (Model  M3A-LP,  MPR  0.06) 
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Figure  68  Transverse  Stress,  v'v'/uj  (Model  M3A-LP,  MPR  0.06) 
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Figure  71  Transverse  Stress,  v'v' / uj  (Model  M3A-LP,  MPR  0.03) 
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Figure  73  TKE,  (7~2 +V1  )/2uJ  (Model  M3A-LP,  MPR  0.06) 
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